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Abstract 

In 1995, over six hundred thousand licensed pilots flew nearly thirty-five million flights 
into over eighteen thousand U.S. airports, logging more than 519 billion passenger miles. 
Since demand for air travel has increased by more than 50% in the last decade while ca- 
pacity has stagnated, congestion is a problem of undeniable practical significance. In this 
thesis, we will develop optimization techniques that reduce the impact of congestion on 
the national airspace. We start by determining the optimal release times for flights into 
the airspace and the optimal speed adjustment while airborne taking into account the 
capacitated airspace. This is called the Air Traffic Flow Management Problem ( TFMP ). 
We address the complexity, showing that it is NP-hard. We build an integer progr amm ing 
formulation that is quite strong as some of the proposed inequalities are facet defining for 
the convex hull of solutions. For practical problems, the solutions of the LP relaxation 
of the TFMP are very often integral. In essence, we reduce the problem to efficiently 
solving large scale linear programming problems. Thus, the computation times are rea- 
sonably small for large scale, practical problems involving thousands of flights. Next, we 
address the problem of determining how to reroute aircraft in the airspace system when 
faced with dynamically changing weather conditions. This is called the Air Traffic Flow 
Management Rerouting Problem (TFMRP) .We present an integrated mathematical pro- 
gramming approach for the TFMRP, which utilizes several methodologies, in order to 
minimize delay costs. In order to address the high dimensionality, we present an aggre- 
gate model, in which we formulate the TFMRP as a multicommodity, integer, dynamic 
network flow problem with certain side constraints. Using Lagrangian relaxation, we 
generate aggregate flows that are decomposed into a collection of flight paths using a 
randomized rounding heuristic. This collection of paths is used in a packing integer pro- 
gramming formulation, the solution of which generates feasible and near-optimal routes 
for individual flights. The algorithm, termed the Lagrangian Generation Algorithm, is 
used to solve practical problems in the southwestern portion of United States in which 
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the solutions are within 1% of the corresponding lower bounds. 
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Chapter 1 


Introduction 


The thesis is structured as follows. In Chapter 1, we introduce two dynamic flow man- 
agement problems that are the topic of this thesis and review the literature to provide a 
framework for our contribution. In Chapter 2, we formally introduce the air traffic flow 
management problem, present our formulation, address the complexity, discuss some 
modeling variations and examine the theoretical properties of our formulation. We prove 
that the proposed constraints are facet defining which provides insight into the computa- 
tional performance. In Chapter 3, we report computational results for the air traffic flow 
management problem and an important special case called the multiple airport ground 
holding problem. In Chapter 4, we discuss the air traffic flow management rerouting 
problem and present the Lagrangian Generation Algorithm. In Chapter 5, we present 
computational results for the air traffic flow management rerouting problem based on 
real data. In Chapter 6, we include some concluding remarks and directions for future 
research. Technical proofs are provided in the appendices. 

1.1 The Air Transportation Industry 

In 1995, over six hundred thousand licensed pilots flew nearly thirty-five million flights 
into over eighteen thousand U.S. airports, logging more than 519 billion passenger miles. 
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Since demand for air travel is increasing by more than 50% each year and capacity is 
stagnating, congestion is a problem of undeniable practical significance. For U.S. airlines, 
the expected yearly cost of the resulting delays was estimated at $3 billion in 1995. In 
order to put this number in perspective, the total reported losses of all U.S. airlines 
amounted to approximately $2 billion in 1991 and $2.5 billion in 1990. In fact, it was 
not until 1995 that the total net profit actually became positive. Furthermore, every day 
700 to 1100 flights are delayed by 15 minutes or more. European airlines are in a similar 
plight. 

Faced with the realities of congestion, the FAA has been using ground-holding policies 
to reduce delay costs. These short-term policies consider airport capacities and flight 
schedules as fixed for a given time period, and adjust the flow of aircraft on a real-time 
basis by imposing “ground holds” on certain flights. Such a flight is then held on the 
ground at its departure airport even if it is otherwise ready for takeoff. 

The basis for ground-holding relies on the fact that while a flight is airborne it incurs 
fuel, safety and other costs that are not applicable before the flight takes off. Thus, 
airborne delays are much costlier than ground delays. If an aircraft departs on time, 
only to encounter airborne congestion as it awaits landing clearance at the destination 
airport, it may incur an airborne delay. However, by delaying its departure with a ground 
delay, the aircraft arrives at its destination at a later time when minimal congestion is 
expected, thus, avoiding the costly airborne delay. Therefore, the objective of ground- 
holding policies is to “translate” costlier anticipated airborne delays to the ground. 

Unfortunately, merely imposing ground holds only addresses a small portion of the 
problem of alleviating congestion. Greater benefit could be realized by taking a more 
complete view of the national airspace. We can control the time that a flight reaches the 
arrival airport by controlling its speed throughout its route, rather than forcing the flight 
to wait on the ground before departure. This may be an attractive option especially if 
the departure capacity is decreasing with time. Moreover, we could adjust the route of 
a flight. For instance, if there is a bad weather system that creates an airspace region of 
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limited capacity, we could reroute a flight along a different path to attempt to avoid the 
weather system. 

1.2 Dynamic Flow Management Problems 

In this thesis, we develop optimization techniques that encompass all aspects of the 
national airspace in order to minimize cost associated with delay. In particular, we 
address two important problems in air traffic control, the Air Traffic Flow Management 
Problem (TFMP) and the Air Traffic Flow Management Rerouting Problem (TFMRP). 

Besides determining release times for aircraft (ground-holding), the Air Traffic Flow 
Management Problem (TFMP) also determines the optimal airspeed adjustment of air- 
craft for a network of airports taking into account the capacitated airspace. Thus, the 
TFMP determines how to control a flight throughout its duration, not simply before its 
departure. This thesis makes the following contributions to this problem. 

We build an integer programming formulation whose objective is the minimization 
of delay costs. The formulation is quite strong as some of the proposed inequalities are 
facet defining for the convex hull of solutions. We address the complexity of the TFMP 
and show that it is NP-hard. We illustrate how our models can be adjusted to account 
for several variations in the problem’s characteristics. When modified for a special case 
previously addressed in the literature, called the multiple airport ground holding problem, 
we prove that the LP relaxation bound of our formulation is at least as strong as all others 
proposed in the literature. For practical problems, the solutions of the LP relaxation of 
the TFMP are very often integral, so there is no need to branch and bound. In essence, 
our formulations reduce the problem to efficiently solving large scale linear programming 
problems. As a result, the computation times axe reasonably small for large scale, realistic 
size problems involving thousands of flights. Short computational times and integrality 
properties are particularly important, since these models are intended to be used on-line 
and solved repeatedly during a day. 
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If we add the final complication, rerouting of flights due to drastic fluctuations in 
the available capacity of airspace regions, we obtain the Air Traffic Flow Management 
Rerouting Problem (TFMRP). In this problem, a flight may be rerouted through a dif- 
ferent flight path in order to reach its destination if the current route passes through a 
region that has very low capacity for reasons usually related to poor weather conditions. 
This thesis will make the following contributions to this problem. 

We present an integrated mathematical programming approach for the TFMRP, 
which utilizes several methodologies, for the problem of minimizing delay costs. In order 
to address the high dimensionality, we begin by presenting an aggregate model, in which 
the problem is formulated as a dynamic, multicommodity, integer network flow prob- 
lem with certain side constraints. Using Lagrangian relaxation, we generate aggregate 
flows that are decomposed into a collection of flight paths for individual aircraft using a 
randomized rounding heuristic. This collection of paths is then used in a packing inte- 
ger programming formulation, the solution of which generates feasible and near-optimal 
routes for individual flights. The overall algorithm, termed the Lagrangian Generation 
Algorithm, is used to solve real problems in the southwestern portion of United States. 
The solutions returned by our algorithm for practical problems are within 1% of the 
corresponding lower bounds. 

Currently, the FAA implements a national ground-holding policy. This policy consists 
of a two-stage process. First, the FAA determines the arrival slots that each airline will 
receive using a first-come, first-serve rule. Next, each airline decides which of their 
scheduled flights to assign to each arrival slot. From this, the FAA calculates the amount 
of ground holding that each flight will incur. In the last decade, several models have 
been developed that use optimization techniques to improve upon the current practices 
for selecting ground holds. The FAA has investigated the possibility of implementing 
these techniques. In the next section, we will discuss how the TFMP and the TFMRP 
fit into the framework of previously addressed problems. 
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1.3 Literature Review 


In Odoni (1987), the problem of scheduling flights in real time in order to minimize 
congestion costs was first conceptualized and introduced. Since then several models have 
been proposed for solving different versions of this problem. The first and simplest version 
considers a single airport and makes decisions about the ground-holds for this Single- 
Airport Problem (SAGHP ). The Multi- Airport Ground- Holding Problem (MAGHP ) was 
the next problem to be addressed. It makes ground-holding decisions for an entire network 
of airports. Thus, the SAGHP and the MAGHP are distinguished by whether delays are 
assumed to propagate in the network of airports as aircraft perform consecutive flights. 

As discussed above, the Air Traffic Flow Management Problem (TFMP) further de- 
termines the optimal airspeed adjustment of aircraft for a network of airports taking 
into account the capacitated airspace. For the Air Traffic Flow Management Rerouting 
Problem (TFMRP), a flight may be rerouted through a different flight path in order to 
reach its destination if the current route passes through a region that has reduced ca- 
pacity primarily due to poor weather conditions. In order to describe the work on these 
problems we consider the following modeling variations: 

1. Deterministic versus stochastic models, which are distinguished by whether the 
capacities of the system (airports and sectors in the airspace) are assumed deter- 
ministic or probabilistic. 

2. Static versus dynamic models, which are distinguished by whether or not the solu- 
tions are updated dynamically during the day. 

The deterministic SAGHP (both static and dynamic) was first formulated as a net- 
work flow problem in Terrab and Odoni (1991). The stochastic SAGHP was formulated 
and solved as a stochastic programming problem in Richetta and Odoni (1993) (static 
case) and Richetta and Odoni (1994) (dynamic case). A review of optimization models 
for the SAGHP is given in Andreatta, Odoni and Richetta (1993). 
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The deterministic MAGHP was formulated as a 0-1 integer programming problem 
in Vranas, Bertsimas and Odoni (1994a) (static case) and in Vranas, Bertsimas and 
Odoni (1994b) (dynamic case). In this thesis and in Andreatta and Tidona (1994), new 
formulations for the MAGHP are proposed. These three models are computationally 
compared in Andreatta and Brunetta (1995) which concludes that our model performs 
the best computationally. We will further prove in this thesis that our model is the 
strongest formulation of MAGHP. A heuristic approach to solving the MAGHP is given 
in Brunetta, Guastalla, and Navazio (1996). Terrab and Paulose (1993) address the 
stochastic MAGHP as a stochastic programming problem. The papers by Matos and 
Ormerod (1995) and by Vranas (1995) discuss a problem in the European network. These 
papers address the “slot allocation” problem. Vranas (1995) shows that this is equivalent 
to the MAGHP in the case where congestion may only arise at the destination airport. 
Ball (1993) and Milner (1995) address the problem of banking flights and Gilbo (1993) 
addresses the problem of dependent arrival and departure runways. 

In this thesis, we present a 0-1 integer programming model for the deterministic, 
multi-airport TFMP which addresses capacity restrictions on the en route airspace. Si- 
multaneously with our work, models addressing enroute capacities were also introduced 
by Lindsay, Boyd and Burlingame (1993). They propose integer programming formula- 
tions for a version of TFMP that tracks a flight as it passes from fix to fix in the airspace. 
The fixes are points in the airspace, not regions of airspace. So the flights only experience 
capacity restrictions at airports and at fixes. As the linear programming relaxations of 
these formulations are not very strong, branch and bound is needed to generate integral 
solutions. However, by developing a wide array of novel formulation strengthening tech- 
niques, the dependence on “pure” branch and bound, as well as the computation times, 
are actually reduced. Helme (1994) has presented a method for the TFMP by designing 
a multicommodity minimum cost flow model over a network in space- time. This method 
has not as yet been fully tested, but it is expected that there will be severe dimensionality 
problems. 
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The problem of dynamically rerouting aircraft has not been addressed to the best 
of our knowledge in the literature. The backbone of our approach is the dynamic net- 
work flow formulation. Ford and Fulkerson (1958) first introduced a dynamic maximum 
flow problem as a standard network generalized to include traversal times between nodes. 
There are algorithms for problems that work with the dynamic network directly. Ford and 
Fulkerson (1958) present an algorithm that solves the dynamic maximum flow problem. 
Wilkinson (1971) and Minieka (1973, 1974) both present an exponential time algorithm 
to solve the universally maximum dynamic flow problem. In addition to these early pa- 
pers, there has recently been considerable research activity on theoretical approaches to 
dynamic network flow problems. Hoppe and Tardos (1994, 1995a, 1995b) have described 
several polynomial time algorithms for discrete dynamic network problems including ap- 
proximate universally maximum dynamic flows, lexicographically maximum flows and 
dynamic transshipment. Fleischer and Tardos (1996) furthered these algorithms by look- 
ing at the analogous continuous-time problems. For a complete review of work done 
on dynamic network flows see the survey papers of Aronson (1989), Bookbinder and 
Sethi (1980) and Powell. Jaillet and Odoni (1995). These advances are not directly rel- 
evant for our problem as our formulation is both multicommodity, integer and involves 
complicating side constraints. 
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Chapter 2 


The Air Traffic Flow Management 
Problem 


In this chapter, we discuss the air traffic flow management problem. In Section 2.1, we 
present an IP formulation for the traffic flow management problem and examine the size 
of the formulation. In Section 2.2, we make the connection with the ground-holding 
problem, by showing how the TFMP formulation can be reduced to model the MAGHP. 
Furthermore, the resulting MAGHP formulation is compared with others proposed in the 
literature. In Section 2.3, we provide some analysis of the polyhedral structure of the 
linear relaxation of our IP formulation of the TFMP. In Section 2.4, the TFMP is proven 
to be NP-hard. In Section 2.5, we extend the TFMP formulation to take into account 
several variations of the model, which incorporate the dependence between arrival and 
departure capacities, hub connectivity, banks of flights and rerouting of aircraft. 


2.1 The 0-1 IP Formulation 

The national airspace is divided into sectors. A map of the United States that displays 
all of the high level sector boundaries is given in Figure 2-1. Each flight passes through 
contiguous sectors while it is en route to its destination. There is a restriction on the 
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number of airplanes that may fly within a sector at a given time. This number is depen- 
dent on the number of aircraft that an air traffic controller can manage at one time, the 
geographic location and the weather conditions. We will refer to the restrictions on the 
number of aircraft in a given sector at a given time as the en route sector capacities. 

Consider a set of flights, T = {1, . . . , F}, a set of airports, K = {1, . . . , A'}, a set 
of time periods, T ={1,...,T}, and a set of pairs of flights that are continued, C = 
{(/',/) : f is continued by flight /}. We shall refer to any particular time period t as 
the “time t.” The problem input data are given as follows: 

Data: 


N f 

Pf 

Dk(t) 
A k (t ) 
Sj(t) 

df 

r f 

s f 

4 

c / 

h,j 

T f 

T) 

T f 


number of sectors in flight /’ s path, 

f 

the departure airport, if i = 1, 

< the (?' — l) s ‘ sector in flight /’ s path, if 1 < i < Nf, 

the arrival airport, if i = Nf, 

departure capacity of airport k at time t, 

arrival capacity of airport k at time t, 

capacity of sector j at time t , 

scheduled departure time of flight /, 

scheduled arrival time of flight / , 

turnaround time of an airplane after flight /, 

cost of holding flight / on the ground for one unit of time, 

cost of holding flight / in the air for one unit of time, 

number of time units that flight / must spend in sector j, 

set of feasible times for flight / to arrive to sector j = [T J f,T J f\, 

first time period in the set Tj, 

last time period in the set Tj. 
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Note that by “flight” , we mean a flight leg between two airports. Also, flights referred 
to as “continued” are those flights whose aircraft relies on an aircraft that has just 
completed a previous flight. 

Objective: The objective in the TFMP is to decide how much each flight is going to be 
held on the ground and in the air in order to minimize the total delay cost. 


We model the problem as follows: 
Decision variables: 



1 , 

0 , 


if flight / arrives at sector j by time t, 
otherwise. 


Note that the wj t are defined as being 1 if flight / arrives at sector j by time t. This 
definition using by and not at is critical to the understanding of the formulation. Also 
recall that we have also defined for each flight a list P f including the departure airport, 
the pertinent sectors and the arrival airport, so that the variable wj t will only be defined 
for those elements j in the list Pf. Moreover, we have defined Tj as the set of feasible 
times for flight / to arrive to sector j , so that the variable wj t will only be defined for 
those times within Tj. Thus, in the formulation whenever the variable t uj t is used, it 
is assumed that this is a feasible ( f,j,t ) combination. Furthermore, one variable per 
flight-sector pair can be eliminated from the formulation by setting w3 ff f = ^ Since 
flight / has to arrive at sector j by the last possible time in its time window, we can 
simply set it equal to one as a parameter before solving the problem. 

To ensure the clarity of the model, consider the following example pictured in Figure 
2-2 which depicts two flights traversing a set of sectors. In this example, there are two 
flights, 1 and 2, each with the following associated data: 


P 1 = (1,A,C,D, E, 4) and P 2 = (2, F, E, D, B, 3). 
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T 

{ 


Figure 2-2: Two possible flight routes. 

If we consider the current position of the aircraft at time f, indicated by the position of 
the airplane icons in Figure 2-2, then the decision variables for these flights at this time 
are given by: 

= l,iu£ t = l,w° t = = 0, and 

w h = = i = o,«* = o, w| it = o. 

Having defined the variables wj t we can express several quantities of interest as linear 
functions of these variables as follows: 

1. The variable u 1 * 3 f t = 1 if flight / arrives at sector j at time t and 0 otherwise, can 
be expressed as follows: 

u 3 f t = w 3 f t - w 3 f t _ 1 and vice versa, ^ u 3 fv . (2.1) 

As expressed earlier, the variables w j f t are only defined in the time range Tj, so 
that = 0. Furthermore, the constraint that a flight must arrive at sector 

j at some time t, originally expressed by the restriction u /,« — 1 can now 

be replaced by the simpler expression w 3 -= = 1. As previously mentioned, this 
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can be handled as a parameter before the problem is solved, thus eliminating many 
variables and constraints. This substitution is fundamental to the performance of 
this model. 

2. Noticing that the first sector for every flight represents the departing airport, the 
total number of time units that flight / is held on the ground can be expressed as 
the actual departure time minus the scheduled departure time, i.e., 

9f= Y ~ d f- 

{ter}, {ter;, 

k=P(f, i)> fc=P(/,D} 

3. Noticing that the last sector for every flight represents the destination airport, the 
total number of time units that flight / is held in the air can be expressed as the 
actual arrival time minus the scheduled arrival time minus the amount of time that 
the flight has been held on the ground, i.e., 

a/= Y tu ),t- r I ~ 9} = Y t ( wk f,t- w f,t-i)~ r f-9f- 

{ter*, (ter*, 

k=P(f,N/)} k=P(f,N f )} 

The objective function 

The objective of the formulation is to minimize total delay cost. Using the variables 
< 7 / and a/ for the amounts of ground and air delay respectively, as defined in items 2. 
and 3. above, the objective function can be expressed simply as follows: 

Min Y Wf9f + C °f a f ]• 

{/e^} 

Substituting the expressions we derived in items 2 and 3 above for the variables wj t 
we obtain the following expression: 
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(TFMP) 

IZ TFMP = Min £ (cf-cj) 
{/6^} 


c / £ 
{<67*. 

*«#»(/,*,)} 

subject to 


{/:P(/,l)=fc} 


£ (w/,t - < Mt), 

{f-.P(f,N f )=k} 


£ (w/,t - w it) < ^(t), 

P(/, *+!)=/, «AT/} 


w /!t+b,, ~ w /,« - °» 


£ <(«/,* - ™), t - 1) + 
(ter;, 

*=*>(/, i)} 


+(c/ - cf )rf/ - cjr/] 


VkelC,teT, 

(2.2) 

vfceK,<€T, 

(2.3) 

Vj 6 6 T, 

(2.4) 

f vfeJ 7 ,teTj,j=P(f,i), 
\ j'=P(f,i+ l),i < Nf, 

(2.5) 
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( 2 . 6 ) 




< 0 , 


' v(/',/) eC,t e Tf, 

\ k = P(f, 1) = P(f, Nf), 


w f, t ~ w f,t- 1 - v/ € T,j G Pf,t eTj, 


(2.7) 


w },t e V/G ^ ,j G Pf,t G Tj. 

The first three constraints take into account the capacities of various aspects of the 
system. The first constraint ensures that the number of flights which may take off from 
airport k at time t , will not exceed the departure capacity of airport k at time t. Likewise, 
the second constraint ensures that the number of flights which may arrive at airport k 
at time t, will not exceed the arrival capacity of airport k at time t. In each case, the 
difference will be equal to one only when the first term is one and the second term is 
zero. Thus, the differences capture the time at which a flight uses a given airport. The 
third constraint ensures that the sum of all flights which may feasibly be in sector j at 
time t will not exceed the capacity of sector j at time t. This difference gives the flights 
which are in sector j at time t , since the first term will be 1 if flight / has arrived in 
sector j by time t and the second term will be 1 if flight / has arrived at the next sector 
by time t. So, the only flights which will contribute a value of 1 to this sum are those 
flights that have arrived at j and have not yet departed from j by time t. 

Constraints (2.5) represent connectivity between sectors. They stipulate that if a 
flight arrives at sector f by time t + If j, then it must have arrived at sector j by time t 
where j and f are contiguous sectors in flight /’ s path. In other words, a flight cannot 
enter the next sector on its path until it has spent l fj time units (the minimum possible) 
traveling through sector j, the current sector in its path. 

Constraints (2.6) represent connectivity between airports. They handle the cases in 
which a flight is continued, i.e. , the flight’s aircraft is scheduled to perform a later flight 
within some time interval. We will call the first flight /' and the following flight /. 
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Constraints (2.6) state that if flight / departs from airport k by time t, then flight /' 
must have arrived at airport k by time t — Sf>. The turnaround time, Sf, takes into 
account the time that is needed to clean, refuel, unload and load, and further prepare 
the aircraft for the next flight. In other words, flight / cannot depart from airport k, 
until flight /' has arrived and spent at least Sf time units at airport k. 

Constraints (2.7) represent connectivity in time. Thus, if a flight has arrived by time 
t, then w J f t , has to have a value of 1 for all later time periods, t' > t. 

Important Remark 

The major reason we used the variables v as opposed to the variables u 3 ft is that the 
former variables nicely capture the three types of connectivity in TFMP: connectivity 
between sectors, connectivity between airports and connectivity in time. Of course, 
given that the two sets of variables are linearly related, the same constraints can be 
captured using the u 3 ft variables. We feel, however, that the variables w 3 f l not only take 
connectivity naturally into account, but also they define connectivity constraints that 
are facets of the convex hull of solutions (see Section 2.3). As we report in Section 3.2, 
the LP relaxation of (TFMP) is almost always integral, i.e., the given formulation is 
a particularly strong one. We believe that the key for this is the use of the decision 
variables w 3 j t in the formulation. 

Size of the Formulation 

Let 


D = max IT*' I, 

{ferjePf}' 3 

be the maximum cardinality of the set of feasible times for flight / to be in sector j taken 
over all / and j, and let 


X = max Nf 

{/€JF} 

be the maximum number of sectors that a flight passes through along its route, taken 
over all flights. Note that X > 2, since the departure and arrival airports are always 
counted as sectors on a flight’s path. Let \T\ be the total number of flights, |X| be 
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the total number of time periods, \K\ be the total number of airports, \J\ be the total 
number of sectors, and \C\ be the total number of flights that are continued. 

The actual number of variables wj t is 


£ £ \t’,\ 

f IKP,) 

since each flight has a different number of sectors and number of feasible time intervals 
associated with it. An upper bound on the number of variables t uj t will be 

\T\ D X. 


The exact number of constraints is 

2\)C\\T\ + \J\\T\ + 2 £ £ 1171 + £ min{|T;i,|r;,|} 

{/e^} {] €P f } {(/',/)ec, 

«=W), 

a =P(f\Nf,)} 

An upper bound on the number of constraints can then be calculated as 

2|/C||T| + \J\\T\ + 2\T\DX + \C\D. 

In order to get a feeling of the size of the formulation let us consider an example that 
adequately represents the U.S. network: 

1. K = 20 representing the most congested airports in the U.S. 

2. |T| = 14 * 12 = 168, representing a 14 hour day with 5 minute intervals. 

3. \J\ = 200, representing 200 sectors. 

4. \T\ = 10000, representing approximately half of the number of daily flights of major 
carriers. 
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5. \C\ = 8000, representing an 80% connectivity among flights. 

6. D = 6, representing an upper bound of half an hour that a flight can be late to any 
given sector. 

7. X = 5, representing an upper bound of at most 5 sectors in a flight’s path. 

For this example the number of variables is at most 300,000 and the number of 
constraints is at most 688, 320. The critical quantities that significantly affect the number 
of variables and constraints are D, X, and \T\. If for example any of these parameters 
doubles, the number of variables doubles and the number of constraints nearly doubles. 

2.2 The Multi- Airport Ground-Holding Problem as 
a Special Case 

As mentioned in Chapter 1, the ground-holding problem is a special case of the TFMP. If 
we remove the sector capacity constraints and the variables associated with the sectors, 
we obtain a new formulation of the MAGHP which, as we demonstrate in Chapter 3.2, 
leads to significant computational advantages compared to alternative formulations that 
have previously been proposed (see Section 1.3). Notice that N f = 2 for all f € X, since 
a flight’s path consists solely of the departure and arrival airports. 

Let us redefine the variables as: 

2// >t = Wjj, for the departure airport, k = P(f, 1), 

Zfj = w k j t , for the arrival airport, k = P(f, 2). 

Also, let Tj be the set of feasible departure times for flight / and let T“ be the set of 
feasible arrival times for flight /. 

Using the new variables, the formulation (TFMP) specializes to the following new 
formulation of (MAGHP): 

(MAGHP) 
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IZmaghp — Min ^ 1 (c/ - C/) 

{/eF} 

E %/.« - Vf.t-i) 

per/} 


+ c / t(z u -z f ,t-i) 

{ter;} 

+{c a f - Cj)df - c a f r f _ 


subject to (yf,t ~ Vf,t- 1) < Dk{t), 

{fterf} 

Vk g /C,t g T, 

(2.8) 

E (z u -z Lt ^) <A k (t), 

V/c G /C, t G T, 

(2.9) 

z f,t ~ yf,t-(r f -df) 5: 0, 

v/ef,i€T;, 

(2.10) 

yf,t - Zf,t—8fi < 0, 

V(/',/)GC,<GT/, 

(2.11) 

1 

7 

IV 

o 

V/ £ iF,t £ Tj, 

(2.12) 

z f,t - > 0, 

v/gt.igt;, 

(2.13) 

y f ,t,Zf,t e {0, 1}, 

V/ G G T. 



The first two constraints incorporate the capacity restrictions of the departure and 
arrival airports. The next constraint is the sector connectivity constraint, which is equiv- 
alent to constraint (2.5) in the TFMP formulation. However, for the ground holding 
problem the only elements in the path are the departure airport and the arrival airport. 
So this constraint connects these two elements by making sure that flight / cannot arrive 
at time t unless it has departed by at least t minus the minimum flight time. The next 
constraint is the flight connectivity constraint, which is equivalent to constraint (2.6) 
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in the TFMP formulation. The last two constraints are time connectivity constraints, 
which are equivalent to constraint (2.7) in the formulation (TFMP). 

Using the previous definitions, an upper bound on the number of variables is 2\F\D 
and an upper bound on the number of constraints is 2|/C||T| + 2>\F\D + \C\D. For the 
same example as in the end of the previous subsection, an upper bound on the number 
of variables in the above formulation is 120,000 and an upper bound on the number of 
constraints is 234, 720. 

If we remove the constraint (2.11) and consider the set /C to be the singleton set, then 
we have a valid formulation for SAGHP, which we will call (SAG HP). We define the 
feasible regions for the formulations (TFMP), ( MAGHP ), and (SAGHP) as IPtfmp, 
IPmaghp and IPsaghp respectively. 

The variables used in the formulation in Vranas et al. (1994a) are defined differently: 
u /, t = 1 if flight / takes off at time t and ?; /)t = 1 if flight / arrives at time t. These are 
linearly related to variables yj t and z^ t as per the relationship given by (2.1). As already 
mentioned, the ground-holding delays can be expressed in terms of these variables in the 
following manner: 

9f=Z tu fA ~ d f , 

{ter/} 

as can the airholding delay, 


a f = tv f,t ~ r f~ 9} ■ 

{ter/} 

In Vranas et al. (1994a), it is assumed that when the departure capacity is large, without 
loss of generality, a f = 0, thus implying that all of the delay would be taken on the ground 
before departure. This gives an equivalent expression for gf as, g/ = 52 teT a tvj >t — ry, 
which contains no departure information, thus eliminating the variables uj tt from the 
formulation. Moreover, instead of the flight connectivity constraints (2.11), the following 
constraints, 

9f ~ (df — Sfi — r f>) < g f , (2.14) 
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establish connectivity between the arriving flight /' and the departing flight / by forcing 
the amount of ground-hold for flight / to be at least the amount that flight /' arrives 
late, g f >, minus the amount of slack time, df - Sf - r f >. The description of the feasible 
space in Vranas et al. (1994a) expressed in the z />t space as per the relationship (2.1) is 
as follows: 


IPvbo = \ z f,t € {O’ 1} y(z/,t ~ ) < Ak{t), y ( z f‘ l z fd-\) — 

9f=J2 ~ z f,t- 1 ) - r f ,g r - {d f - s r - r r ) < g f ,z u - > oj. 

{ter;} 


Terrab and Paulose (1993) use the same variables, v ft as in Vranas et al. (1994a). 
However, they express the flight connectivity constraints as follows: 


y vu — y v f,v < o. ( 2 .i 5 ) 

{ter f a ,t<T} {t’zrp, 

t'<T—s ft —(rj—df)} 

Constraint (2.15) forces connectivity, since if the second sum is zero then flight /' has 
not landed by time r — sj> — ( r f — df ), which is time period r minus the turnaround 
time, minus the flight time of /. This forces the first sum to be zero so that flight / can 
not land before time r. The description of their formulation expressed in the Z/ if space 
as per the relationship (2.1) is: 


IP T p = [zf t € {0,1} YXz t ,t ~ z f,t-i) < Ak{t), 5Z ( z /,t 2 /,t- 0 ~ 

1 / {ter/} 


y z f,t-\) — y ( z f,t' — Zf\t>~i) < o, Zf^ — Zf it ~i > oi. 

{tery.Kr} } 

t’<T-Sft-{Tf-df)} 


If we specialize our formulation for the case of large departure capacities and use only 
the variables, Zf it (j//,t = Zf,t+( rf -d f )), we obtain: 
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IP'maghp — \ z f,t £ { 0 > 1 } ( z f t ~ z U-i) — Ak(t), ^2 (z u — Zf <t - 1 ) = 1 , 

1 {/:ter;} {ter;} 

Z f,t+( r f—df) z f,t z f,t— 1 — O}. 

In all of these formulations, the expression ^2 t £Tf( z f,t ~ z f,t- 1 ) = 1 reduces to the 
expression Zjj f = 1. This telescoping property is due to the unique definition of the 
decision variables as flights arriving by some time t rather than at time t. 

If we denote the polyhedra corresponding to the linear programming relaxations of 
IP'maghpi IPvbo, and IPtp as P'maghpi Pv bo, and P TP and denote their corresponding 
values as Z' MAGH p, Zvbo > and Z then we can state the following proposition whose 
proof is included in Appendix A. 

Proposition 1 IP TP = IP VB o = IP'maghp ^ P'maghp ^ p tp Q Pvbo and corre- 
spondingly , Zy bo — Z T p ^ ZftjAGHP — ^ Z MAGHP ~ IZyBO = IZfp. 

Therefore, the LP relaxation of {MAGHP) gives bounds that are at least as strong 
as those from the LP relaxations of either Vranas et al. (1994a) or Terrab and Paulose 
(1993). 

2.3 Insights from the Polyhedral Structure 

In Section 3.2, we report computational results for the TFMP based on the formulation 
( TFMP ). Even for large scale problems and for a variety of problem parameters the 
solutions of the LP relaxation of both {TFMP) and {MAGHP) were integral. In the 
tradition of polyhedral combinatorics in mathematical programming, we examine the 
polyhedral structure of Ptfmp a nd Pmaghp in order to obtain a deeper understanding 
of why this formulation performs so well computationally. Given a set S we denote 
with conv{S) the convex hull of solutions in S. In particular, we address the following 
questions: 
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1. Are the polyhedra P T fmp and Pmaghp integral? If not, is the optimal solution to 
the optimization problem integral if we impose the simplification that c g = and 
c a = cj for all / € T7 

2. Are the constraints in ( TFMP ) and ( MAGHP ) facets of con.v(I P TFM p) and 
conv(I Pmaghp) respectively? 

We summarize our findings in the following theorem: 

Theorem 2 

1. The polyhedra Ptfmp and Pmaghp are not integral. Even with the simplification 
that c g = cj and c a = c) for all f € T , integral solutions are not obtained. 

2. Inequalities (11), (12), (13) and (14) are facet defining for cov.v(IP M aghp), while 
the constraints (9) and (10) need not be. Inequalities (5), (6) and (7) are facet 
defining for conv(I Ptfmp)- while the constraints (2), (3) and (4) need not be. 

As the proofs of the theorem are somewhat technical, we have placed them in appen- 
dices B, and C respectively. 

The previous theorem gives some partial insight on the usefulness of the new variables 
we introduced, which make it easy to express sharply the various types of connectivity 
in the problem. While the formulations are not integral, the inequalities that the three 
types of connectivity impose are indeed facet defining. As the solutions obtained were 
integral for a wide spectrum of examples and parameters, we did not investigate further 
the determination of other facets. 

2.4 Complexity of the TFMP 

In this section we show that the TFMP is an NP-hard problem. 

Theorem 3 The TFMP with all capacities equal to 1 is NP-hard. 
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Proof. We show that job-shop scheduling (see Garey and Johnson (1979)) reduces to 
TFMP. 

JOB SHOP SCHEDULING PROBLEM (JSP) 

INSTANCE: Number m G Z + of processors, set J of jobs, each j G J consisting of an 
ordered collection of tasks t h [j], 1 < k < rij, for each task t a length l(t) G Zq and a 
processor p(t) G {1,2, . . . ,m}, where p{t k [j}) ^ p{t k +i[j]) for all j € J and 1 < k < n h 
and a deadline D G Z + . 

QUESTION: Is there a job-shop schedule for J that meets the overall deadline, i.e., a 
collection of one-processor schedules u l mapping {t : p{t) = i} into Zq, 1 < i < m, such 
that 


(Ti(t) > <7 x (t') implies a^t) > a^t') -I- l(t), 

Mt k+ i[j\) > <ri(tk\j]) + l(t k [j ]), where i' = p{t k+1 [j]), i = p(t k [j]) 

for all j G J,k. G {1,. 

and rr t (t nj [j]) + l(t nj [j}) < D, where i = p(t Uj [;']) for all j G J. 

For each job we create an aircraft. For each processor we associate an airport or 
sector. Task t k [j) of job j corresponds to a flight segment, f k [j) of aircraft j. Given a 
collection of tasks, t k \j] of job j, we associate a list of airports and sectors to be visited 
by aircraft j. Furthermore, the processing time of task t k [j] corresponds to the time 
required to perform the flight segment, f k [j\. 

We obtain a list of airports and sectors, 


M 1 9 2 A k 

! pj I • • ' > 1 



and a list of the flight segment times, 



i k 

■ - l Ai 
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for each aircraft j by the relationships: 

= 'Wi]). 
t% = '(',(21). 
t% = '(',[3]), 

Aj’ = p(i,[n, 1), 'J, = '(',[ n ,l)- 

So by finding a job-shop schedule that satisfies the given conditions, we will find a 
solution to the transformed problem such that all flights are performed by the deadline 
D. Also, according to the relationship 

Vi'(tk+i[j]) > <Ti(tk\j]) + l(tk[j}) where i' = p(t k+1 [j}),i = p(t k [j]), 

no two tasks will ever performed simultaneously on the same processor, which is equiva- 
lent to limiting the capacities of airports and sectors to one. Moreover, the relationship, 

cr ,(t) > fTj(t') implies rr^t ) > <7i{t') + l(t), 

dictates that a task can not be processed unless the previous task has completed. 

This stipulation guarantees connectivity between flights, and sectors, as specified by 
the set of tasks for each aircraft. Thus, all of the constraints of the TFMP will be satisfied 
if and only if there exists a feasible job-shop schedule. ■ 

2.5 Modeling Variations 

Our goal in this section is to demonstrate that the formulation (TFMP) can be easily 
extended in many directions to take into account several variations of the model. 


A) =p(tj[l]), 

S] ~p(tj[ 2 D, 
Sj =p(tj[ 3 ]), 



2.5.1 Dependence Between Arrival and Departure Capacities 

The interdependence between the arrival and departure capacities of airports results from 
the fact that the same runways are often used for both arrivals and departures. Thus, 
the runway allocation will determine how an airport’s available capacity is allocated 
between the arrivals and departures at a given time. By operating imder a specific runway 
configuration, arrival and departure capacities can be adjusted. This will significantly 
influence airport efficiency. 


22R 



Figure 2-3: Complete Runway Configuration for Logan Airport. 

By choosing a particular configuration of runways for a given time, the capacity 
allocation will be fixed. The complete set of runways for Logan Airport is given in 
Figure 2-3. A common configuration used at Logan Airport is to use runways 4L and 
4R for arriving flights and use runways 9 and 4R for departing flights. Notice that since 
runway 4R is the longest runway and certain types of aircraft require a long runway, it 
is used for both arrivals and departures. 

In Figure 2-4, we have shown a hypothetical runway allocation. Since it takes longer 
for an aircraft to arrive than to depart, if the airport is operating under the allocation 
shown in Figure 2-4 and all the capacity at Logan Airport is allocated to arrivals then 52 
flights could arrive and if all the capacity is allocated to departures then 62 flights could 
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depart within an hour. If only less than 50 flights are allowed the depart within an hour, 
then the capacity allocation is given by the equation, 3 A k (t) + D k (t ) = 3 * 52. Whereas, 
if more than 50 flights are allowed to depart within an hour, the capacity allocation is 
given by the equation, A k (t) + 3 D k (t ) = 3 * 62. 

We review briefly ideas introduced in Gilbo (1993) and Vranas et. al. (1994a). We 
represent the runway allocation by a set of linear constraints indexed by i for airport k 
at time t of the type 


Q kt Dk(t) + (3 l kt A k (t) < Y kt , Vfc e /C,t £ T ,i G I kt , (216) 

where a l kt . /3 l kt , and % t are given constants. In the example of Figure 2-4, there were two 
linear constraints. The set I kt will determine the number of linear constraints for each 
airport k at time t. 


A k «) 



Figure 2-4: Rimway Departure/ Arrival Allocation for a Specific Configuration. 

The region formed by the above constraints gives a complete depiction of all the 
possible runway allocations at a given time, and likewise, all possible departure and 
arrival capacity assignments. 

In order to solve this variation, we treat D k (t) and A k (t) as variables that satisfy 
constraints (2.16) and add them to ( TFMP ). We can further reduce the size of the 
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resulting formulation by eliminating the variables D k (t) and A k (t) by incorporating con- 
straints (2.2) and (2.3) taken at equality into (2.16) as follows: 

4t E ( W U ~ w f,t- 1) + E ( w f,t ~ ™u- 1) ^ 7a- 

{/:te7*fc=P(/, 1 )} {/:ter* 1 fc=P(/,7V / )} 

The addition of this constraint to ( TFMP ) incorporates the dependence between the 
arrival and departure capacity assignments without the addition of any new variables. 

2.5.2 Hub Connectivity with Multiple Connections 

Given that many airlines now control key hub airports through which most of their flights 
are directed, it is no longer obvious which aircraft will fly a subsequent flight. At these 
hubs, many airplanes are capable of performing any one of multiple consecutive flights. 
We refer to the issue of assigning aircraft to continuing flights as hub connectivity. This 
can be achieved by extending the model as follows: 

For each arriving flight /' that is continued there is a set of flights R f > that can continue 
flight /'. Introducing the 0 - 1 variables x f ' f , which take on the value 1 if flight /' is 
continued by flight / G Rf and 0 otherwise, we alter constraint (2.6) as follows 

w h - s 1 - x r.h V(f, f)eC,te T), 

k = P(f, l) = P(f',N,,), 


and add the constraint that each continued flight /' has to be assigned to a flight in Rf. 

E */'./ = !• 

{/£*/'} 

2.5.3 Banks of Flights 

With the evolution of the hub and spoke system, airlines have a set of flights (banks) 
that are scheduled to arrive at a hub airport and another set scheduled to depart within a 
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small time window of the arrival bank. Each arriving aircraft will be assigned to perform 
at most one of the departing flights. This situation is similar to hub connectivity, except 
that the objective of hub connectivity is to minimize delay costs, while the objective of 
the airlines who rely heavily on banking flights is to minimize the spread between the 
arrival times of the first and the last flight in the bank. Let B be the set of flights in a 
bank. We define the decision variables 


VB,t = 


1 if the first flight f in B arrives by time t 
0 otherwise 


Z B. t 


I I if the last flight f in B arrives by time t 
0 otherwise 


These definitions require the constraints: 


ys,t - Wf,t — 0, V / e B,t € Tj , k — P(f , Nf) 

z B ,t ~ w k f t <0, V/ G B,t <E Tf, k = P{f, N f ). 


We also need the additional time connectivity constraints for these variables 


yB,t — VB,t- 1 >o, VfeT 

z B,t — z B,t-i ^ o, wer. 


The objective function of minimizing the “spread” in the arrival times for the flights in 
the bank B can be modeled as follows: 


min J2 - 2b, t— l) ~ ~ VB.t-i)- 

{teT} {teT} 


39 



This is equivalent to determining 


min T 

{ter} 


\ 

max t{w k f i -w) t _ l )- min t {w) t - w) 

\k=P(f,N f )} k=P(f,Nj)} 


'With the addition of these new variables, new constraints, and new objective function, 
banking can be incorporated into the formulation. Alternative approaches to handle 
banking constraints are proposed in Ball (1993) and in Milner (1995). 


2.5.4 Rerouting of Aircraft 

As we already mentioned, very often extreme weather conditions force the capacities of 
some sectors (and airports) in the national airspace system to drop significantly or even 
to become zero. Air traffic controllers are then forced to use alternative routes for aircraft 
passing through these sectors to accommodate these changes in capacities (see Figure 4-1 
for an example). Currently, these rerouting decisions are handled through the experience 
of the air traffic controllers and not through a formal optimization model. Rerouting 
will be discussed extensively in Chapters 4 and 5. We also note that the formulation 
(TFMP) can be extended to accommodate dynamic rerouting decisions. However, these 
extensions quickly lead to very large formulations, especially when considering realistic 
size problems. 

There are two approaches that we may take in attempting to formulate the rerouting 
problem by extending the TFMP formulation. We have called these the path and sector 
approaches. The path approach assumes that each flight has a set of routes to choose 
from in deciding which is optimal. The sector approach makes no assumption about the 
routes, but rather selects that path as a series of possible sectors. 

The path approach first defines Qf as a set of possible routes that flight / may fly. 
In the formulation {TFMP) we have assumed that Qf only contains one route, which 
we have denoted as Pj. In order for the formulation to be of manageable (but still large) 
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size we need to restrict the size of Q f . We extend the TFMP variables in the following 
manner: 


/ 



if flight / arrives at sector j by time t along route r, 
otherwise. 


Clearly, the variables u? f t defined in Section 2.1 can be written as: 



E 


w 


f.f 


{r eQ/} 


Moreover, since the departure and arrival airports will remain the same for a given flight 
over all routes, P(f, 1) and P(f,N f ) will be independent of the particular route. Using 
the newly defined variables we can modify the T FMP to include rerouting with the extra 
restriction that only one path may be used per flight. The size of the resulting formulation 
will be at most a factor max/ \Q } \ larger than the T FMP formulation. This implies that 
we may be able to handle problems with a relatively small number of alternative paths. 

The sector approach decides at each sector in its route which sector to enter next. 
We need to define N(f,j ), the set of sectors that flight / can enter immediately after 
exiting sector j, as well as P{f,j), the set of sectors that flight / can enter immediately 
before entering sector j. We extend the TFMP variables in the following manner: 


{ 1, if flight / arrives at sector j' from sector j by time t, 
0, otherwise. 


Clearly, the variables w^ t defined in Section 2.1 can be written as: 



E 


w 


■ 


0'€N(/,j)} 

As before, the departure and arrival airports will remain the same for a given flight over 
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all routes; P(f, 1) and P(f,Nf ) will be independent of the particular choice of sectors. 
Furthermore, by using the newly defined variables, we can modify the T FM P to include 
rerouting with the necessary constraints that each flight can only travel from one sector 
to one other sector. 
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Chapter 3 


Computational Results for the 
Traffic Flow Management Problem 


In this chapter we report the results of a series of computational experiments that we 
conducted. In Section 3.1, we provide computational results comparing our formulation 
of MAGHP with others in the literature using data from Vranas et al. (1994a). In Section 
3.2, we report computational results of the traffic flow management problem using several 
datasets including some real data provided by the FAA. 

In performing the computational experiments, we aim to address the following ques- 
tions: 

1. How frequently are the solutions obtained by solving the LP relaxations of ( T FMP) 
and {MAGHP) integral? 

2. How is the integrality of solutions affected by the various problem parameters and 
the size of the problem? 

3. How is the computational time required to obtain an optimal solution affected by 
the various problem parameters and the size of the problem? 

4. How does the present approach compare with other approaches in the literature? 
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5. Given that the TFMP needs to be solved on line for controlling air traffic in the US, 
perhaps the most important question to ask is: Can large problems of realistic 
size be solved in reasonable computational times? In other words, is the present 
approach a realistic method to control air traffic in the US? 

3.1 Ground-Holding Problem Computations 

We performed computational experiments on datasets used in Vranas et al. (1994a) on 
the Ground- Holding Problem. Specifically, we looked at the datasets consisting of 2 and 
6 airports with 500 flights per airport, totalling 1000 and 3000 flights respectively. Some 
adjustments in the data were necessary in order to accommodate the differences between 
the two models. In particular, the previous model did not include of any departure data, 
as all of the optimization was done with respect to arrivals. Thus, we generated departure 
data (times and capacities) that were compatible with the existing arrival data. 

As in Vranas et al. (1994a), for each of these cases, four levels of flight connectivity 
were considered. These levels give the ratios of continued flight to total flights, \C\/\T\, 
as 0.20, 0.40, 0.60, and 0.80. 

We considered 15 minute time intervals taken over a 16 hour day. All experiments 
were performed on a Sim SPARCstation 10 model 41. GAMS was used as the modeling 
tool and CPLEXMIP 2.1 was used as the solver. The results that we obtained using the 
above datasets and our (MAG HP) formulation are summarized in Tables 3.1 and 3.2. 


1^1 

\c\m 

Dep Capacity 

Arr Capacity 



1000 

0.20 

32 

15 

262 

0 

1000 

0.40 

17 

10 

741 

4 

1000 

0.60 

20 

14 

359 

0 

1000 

0.80 

20 

20 

283 

mg 


Table 3.1: Results at the infeasibility border for 1000 flights. 
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m 

\c\m 

Dep Capacity 

Arr Capacity 

Time 

% Nonint 

3000 

0.20 

20 

20 

5475 

0 

3000 

0.40 

20 

20 

4703 

0 

3000 

0.60 

20 

20 

5407 

0 

3000 

0.80 

20 

20 

9411 

0 


Table 3.2: Results at the infeasibility border for 3000 flights. 

Tables 3.1 and 3.2 give results at the infeasibility border for each case. The infeasibility 
border is the set of critical values for the departure and arrival capacities, in units of 
flights per time interval, under which the problem becomes infeasible. We expect that it 
is in this region that the problem is very relevant practically and is harder to solve. The 
critical capacity levels were found by a series of trial and error tests. The times reported 
are in CPU seconds and the % Nonint column is the percentage of total flights whose 
solution was noninteger. If we compare these results with the results from Vranas et 
al. (1994a) (see Tables 3.3 and 3.4), we can see that the largest amount of improvement 
occurred in the integrality of the solutions. The computational times for solving our LP 
for 1000 flights (see Table 3.1; column Time) are comparable to the time required to solve 
their LP (see Table 3.3; column LP Time), while for the 3000 flights the LP in Vranas 
et al. (1994a) was solved faster. However, our solutions are for the most part already 
integral (the only instance where the solution was not integral was the 40% connectivity 
instance of the 1000 flight example). The total amount of time required to find an integral 
solution from the LP in Vranas et al. (1994a), found in the total time column, includes 
the time required to solve the LP relaxation, found in the LP Time column, plus the 
time required to perform a Branch & Bound heuristic. If we compare the amount of time 
required to find an integral solution, we see a significant improvement in computational 
time. 

Tables 3.5 and 3.6 were constructed to demonstrate how computational time and in- 
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\C\/\r\ 

Dep 

Capacity 

Arr 

Capacity 

LP 

Time 

Total 

Time 

% Nonint 

1000 

0.20 

oo 

(12,14) 

258 

374 

6.3 

1000 

0.40 

oo 

10 

327 

894 

8.4 

1000 

0.60 

oo 

11 

377 

6958 

12.8 

1000 

0.80 

oo 


lUgkfcM 

9512 

16.8 


Table 3.3: Previous Results at the infeasibility border for 1000 flights. 


m 

ici/m 

Dep 

Capacity 


LP 

Time 

Toted 

Time 

% Nonint 

3000 

0.20 

OO 

12 

1453 

11360 

not given 

3000 

0.40 

oo 

18 




3000 

0.60 

oo 

17 




3000 

0.80 

oo 

18 

3072 

25021 



Table 3.4: Previous Results at the infeasibility border for 3000 flights. 

tegrality are affected by changes in the capacities, i.e., how well does the model perform 
when the capacities are not at the infeasibility border? These results suggest that the 
computational time did not change significantly at different capacity levels. For the one 
case in which the solution was not completely integral, (1000 flights at 40% connectivity), 
increasing the capacity resulted in integral solutions. 
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Dep Capacity 

— 

Arr Capacity 

Obj Value 

Time 

% Nonint, 

1000 

wmM 

32 

17 

50750 

342 

0 

1000 

HU 

32 

16 

55450 

227 

0 

1000 

0.20 

32 

15 

63525 

262 

0 

1000 

0.20 

32 

14 

inf 

- 

- 

lEST 

BEB 

18 

12 

47000 

290 

0 

1 

KH 

18 

10 

79916 

521 

2.2 

1 

0.40 

17 

10 

88241 

741 

4 

tB 

0.40 

16 

10 

inf 

- 

- 


wmm 


18 



0 


1 

| 

15 

33292 


0 



20 

14 

39266 

359 

0 



20 

13 

inf 

- 

- 

1000 


30 

30 


183 

0 

1000 

0.80 

20 

20 


283 

0 

1000 

0.80 

19 

19 

inf 

- 

- 


Table 3.5: Results for varying capacity levels for 1000 flights. 

3.2 Air Traffic Flow Management Problem Compu- 
tations 


We next performed experiments on a connected network of four airports: Boston Logan 
(BOS), NY LaGuardia (LG A), Washington National (DCA) and a node representing 
all other airports (X). Three hypothetical sectors that flights must use before landing 
at LaGuardia, were also introduced into the model. Flights coming from DCA would 
traverse two of these sectors before reaching LGA. Flights coming from BOS would 
traverse one of these sectors before reaching LGA. A certain fraction of the remaining 
flights would use the two sectors approach, while the other flights would enter through 
the one sector. See Figure 3-1. 
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1*1 

|C|/I*I 

Dep Capacity 

Arr Capacity 

Obj Value 

Time 

% Nonint 

3000 

0.20 

30 

30 

42000 

4537 

0 

3000 

0.20 

20 

20 

228000 

5475 

0 

3000 

0.20 

19 

19 

inf 

- 

- 

3000 

0.40 

30 

30 

42000 

5062 

0 

3000 

0.40 

20 

20 

234000 

4703 

0 

3000 

0.40 

19 

19 

inf 

- 

- 

3000 

0.60 

30 

30 

42000 

5629 

0 

3000 

0.60 

20 

20 

234000 

5407 

0 

3000 

0.60 

19 

19 

inf 

- 

- 

3000 

0.80 

30 

30 

42000 

6021 

0 

3000 

0.80 

20 

20 

252000 

9411 

0 

3000 

0.80 

19 

19 

inf 

- 

- 


Table 3.6: Results for varying capacity levels for 3000 flights. 


The three airports (BOS, LGA, DCA) and the three sectors were the only capaci- 
tated elements in the system. The other sectors were allocated unlimited capacity. We 
performed one set of experiments for 200 flights over a 24 hour time period and another 
set for 1000 flights over a 24 hour time period. The 200 flight dataset was obtained from 
the January 1993 Official Airline Guide (OAG). For the larger set of 1000 flights, the 
data was generated by the Pseudo-OAG Generator (POAGG) which is flight generation 
software developed at Draper Laboratories that realistically mimics the flight schedules 
of the OAG. All models were programmed in GAMS, ran on a Sun SPARCstation 10 
model 41 and solved with the solver CPLEXMIP 2.1. For most of the test cases the 
time interval was 5 minutes long. Since some of the sectors could be crossed in under 15 
minutes, we tried to select a time interval that would capture as many sectors as possible 
without becoming prohibitively large. With this in mind, we decided to use a 5 minute 
interval whenever possible. 

For the set of 200 flights, the time frame was 24 hours divided into discrete time 
units of 5 minutes each. To solve the problem CPLEX requires 234 seconds CPU time. 
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Figure 3-1: Sector flow model for test case. 

Moreover, the resulting optimal solution was integral. 

We were able to solve the 1000 flights problem at the infeasibility border over a 24 
hour time period considering 15 minute intervals in 436 seconds CPU time. The optimal 
solution was once again integral. For the complete set of results see Table 3.7. Notice 
that the computation time in CPU seconds varies very little with the capacity restrictions 
in flights per time interval and that the solutions were completely integral. 

Lastly, we obtained two realistic size datasets obtained directly from the OAG flight 
guide. This dataset has also been used to solve similar problems at the MITRE Corpo- 
ration. The first dataset consists of 278 flights, 10 airports and 178 sectors, tested over 
a 7 hour time frame with 5 minute intervals. The second of these datasets consists of 
1002 flights, 18 airports, and 305 sectors tested over an 8 hour time frame with 5 minute 
intervals. 

The sector crossing times, sector and airport capacities, and required turnaround 
times were all provided by the FA A. Nothing used in these datasets was generated or 
hypothesized. We believe that these datasets are very comparable to the problem being 
solved everyday by the FA A. 

For the first problem, consisting of 43226 constraints and 18733 variables, we found 
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Sector Capacity 

Dep Capacity 

Arr Capacity 

Obj Value 

Time 

% Nonint 


20 

20 



0 


20 

20 



0 

20 

10 

10 

inf 

- 

- 

15 

15 

15 

68725 

427 

0 

12 

12 

12 

244225 

450 

0 

11 

11 

11 

inf 

- 

- 

10 

12 

12 

24225 

456 

0 

5 

12 

12 

24350 

432 

0 

4 

12 

12 

250975 

466 

0 

3 

12 

12 

295225 

459 

0 

2 

12 

12 

inf 

- 

- 


Table 3.7: Results for varying capacity levels for 1000 flights. 

an optimal solution in 1141 seconds. Furthermore, the solution obtained was completely 
integral. The second and larger dataset consisting of 151662 constraints and 69497 vari- 
ables, was solved to optimality in 29920 seconds, again achieving completely integral 
solutions. 

In summary, to address the questions we raised in the beginning of this chapter we 
remark: 

1. In all but one instance in MAGHP and all instances of TFMP the relaxations of 
(MAG HP) and (TFMP) were integral. 

2. The integrality of solutions was not affected by problem parameters, nor the size 
of the problem, except for the one instance in which the solution was non-integral. 

3. The computational time required to obtain an optimal solution increases with the 
degree of connectivity as well as with the size of the problem. 

4. Our approach improves upon earlier work particularly in obtaining integral solu- 
tions. This is evident from the fact that our formulation was almost always integral 
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especially in realistic instances. In contrast, for similar test cases (they used fixes 
while we used sectors) the formulation of Lindsay et al. (1993) was not integral, so 
they needed to improve their formulation by using automatic constraint generation 
techniques and also to use branch and bound. 

5. We are able to solve large, realistic size problems in a reasonable amount of time. 
In addition, because we were able to solve the two instances of the TFMP with real 
data, we are very optimistic that our approach can effectively address the TFMP. 
Indeed, the reason we did not solve bigger problems is the difficulty of obtaining 
real data and memory restrictions of the SPARCstation. 
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Chapter 4 


The Air Traffic Flow Management 
Rerouting Problem 


This chapter addresses the traffic flow management rerouting problem (TFMRP), i.e., 
the problem of dynamically rerouting aircraft in the air traffic control system in order to 
avoid airspace regions that have reduced capacities primarily due to bad weather. The 
overall objective being the minimization of delay costs. In Section 4.1, we formulate 
the TFMRP as a integer multicommodity dynamic network flow problem. In Section 
4.2, we discuss the multiple airline problem. In Section 4.3, we present the Lagrangian 
Generation Algorithm for solving this problem. In Section 4.4, we model the TFMP and 
the MAGHP as dynamic network flow problems. 

In the United States, the control of air traffic is centered on 22 regional control 
centers. These centers receive information from aircraft and ground-based radars on 
location, altitude, and speed of aircraft, as well as weather information. When the 
weather conditions are poor, the capacities of some airports and sectors in the national 
airspace are forced to drop significantly or even to become zero. In this situation, the Air 
Traffic Command Center (ATCC) initiates a process to reschedule and reroute flights so 
that the delay costs caused by the weather conditions are kept to a minimum. Aircraft 
must then fly alternative routes if they were scheduled to pass through airspace regions of 
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reduced capacity (see Figure 4-1 for an example). Currently, these rerouting decisions are 
determined through an iterative process between the ATCC and the Airline Operations 
Centers (AOC). The ATCC contacts each airline’s AOC concerning the necessity of 
rerouting. Each AOC then determines a set of new flight paths that they would like to use 
to complete their scheduled flights given the new limited capacity scenario information. 
This problem that each AOC faces when poor weather conditions limit the capacity of 
the national airspace, is the topic of this chapter. 

4.1 The Integer Multicommodity Dynamic Network 
Flow Formulation 

In this section we present an integer, multicommodity dynamic network flow model of 
the TFMRP. There are several components to the model. These include the dynamic 
network, the aggregated flow variables, the non-aggregated variables, and the capacity 
constraints. We will describe each of these in detail below. 

We first describe the dynamic network that models the air traffic system. We create 
a graph which represents the actual geographical picture of the airport/airspace system. 
The nodes of the graph represent both airports as well as sectors. The example in Figure 
4-2 of four airports and six sectors demonstrates how the nodes and arcs of the network 
are constructed. 

The outlined regions are the sectors and the stars are the entrance and exit points 
for the sectors. We define one node for each of these sector crossing points. The circles 
are the airports. Each airport is represented by four nodes as described below. The arcs 
connect the entrance and exit points of a sector as well as these sector crossing points 
and the airports. Each arc (i,j) has a corresponding travel time, t t j . We assume that 
each sector has a limited number of entrance and exit points. In order to represent delay 
in the network we define self-loops with a travel time of one. 

The commodities in the network are defined as origin-destination pairs of airports. 
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Figure 4-2: A network corresponding to four airports and six sectors. 

So, if there are A airports and flights between all airports are flown, then there are exactly 
A(A - 1) commodities. However, if we wanted to distinguish between certain character- 
istics of flights such as aircraft type, we could do this by breaking the commodities down 
even further. We will discuss the multiple airline problem in Section 4.2. 

In order to model airport z we use four nodes i B , zc, and (see Figure 4-3). 



Figure 4-3: An airport is modeled with four nodes. 

Node ia is used to track all the incoming flights to airport z. Once a flight has landed 
at the airport, it is initially at node i A - It can either proceed to node i B or to node i c . 
Node i B represents the situation in which an aircraft has completed all of its required 
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flights for the day. Consequently, the flow into node is is removed from the network for 
the remaining time. Node ic represents the situation in which an arriving aircraft must 
perform at least one more flight during the day. We use k to index the commodities 
of incoming flights and by k' the commodities of outgoing flights. At node i c flow is 
exchanged from commodity k to commodity k' . The delay arc at node ic models the 
situation in which an aircraft arrives before the continued flight is scheduled to depart. 
Flow then proceeds from node i c to node i D . At node i D , new aircraft are introduced 
to the network and all the flights departing from the airport leave from this node. The 
delay arc at this node represents ground holding of flights. 

Let J\f = ( S,S ) be the network formed from airports and sectors, as described above. 
The set of commodities is denoted by {1 ,...,A(A - 1)}, where A is the number of 
airports. We discretize the time horizon into a set of time periods, T = {1, . . . , T}. We 
refer to any particular time period t as the “time t” Note that by “flight”, we mean a 
flight leg between two airports. A flight that is “continued” relies on an aircraft that has 
just completed a previous flight. The problem input data are given as follows. 

Data: 


cm 

r(i) 

orig(k) 

dest(k) 

N(k) 

T 

d f 

c* 


= minimum travel time along arc 
= capacity of sector i at time t , 

= turnaround time required to refuel, reload and clean an aircraft at airport i, 
= airport of origin for a commodity k flight, 

= airport of destination for a commodity k flight, 

= set of arcs that a commodity k flight can use, 

= set of flights, 

= scheduled departure time of flight / 6 T, 

= cost of holding a flight on the ground for one unit of time, 

= average cost of flying an aircraft for one unit of time, 
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H = 
C = 
T } = 
HI) = 
k'(f) = 
Sup k {t ) = 

Dem,k{t) = 


total amount of scheduled flying time for all flights, 
set of flights that are continued, 

set of feasible departure times for flight / where / € C , 
commodity of flight / where / € T , 

commodity of the flight which preceeds flight / where f eC , 
number of flights of commodity k that are scheduled to depart at time t, 
that are not continued, 

number of flights of commodity k that are scheduled to land at time t, 
at the latest, and do not continue to a later flight this day. 


In order to reduce the dimensionality of the problem, we aggregate some of the vari- 
ables over flights. By doing this, we lose the ability to distinguish some of the data 
over the flights. In particular, for the TFMP in Section 2.1, the costs of ground and 
air delay depended upon the flight, cj and c a f . Since we are aggregating over flights in 
this formulation, this is no longer the case. The delay cost variables given above are 
simply c 9 and c a which contain no flight information. Moreover, the turnaround times 
for the TFMP also depended upon the flight, s/, whereas the turnaround times in this 
aggregated formulation can no longer depend upon the flight. Instead, the turnaround 
time depends upon the airport, r(i). 

The most important data that is affected by the aggregation concern the continued 
flights. For the TFMP, each continued flight is assigned a unique flight to precede it. 
In the following formulation, each continued flight is assigned a unique commodity to 
precede it. Thus, for each continued flight, / € C, we define the commodity of flight / as 
k(f) and we define the commodity of the corresponding flight that precedes flight / £ C 
as k'(f). So, instead of forcing flight / 6 C to be a continuation of flight /', we ensure 
that for every continued flight there must be aircraft of commodity k'(f) available for 
flight / € C to use. 


57 



The aggregated variables are defined as: 


Xij(t) = number of flights of commodity k that depart from 
node i at time t and arrive at node j at time t + t hJ . 


Note that these variables are flow variables and not flight variables. So, in order 
to recommend flight paths we must have a method for disaggregating, i.e. for convert- 
ing these flow variables to flight variables. In Section 4.3.2, we propose a method for 
performing the disaggregation, namely, the randomized rounding heuristic. 

We also introduce flight variables for continued flights as follows: 



if the aircraft performing flight / € C is ready for departure at time t, 
otherwise. 


By defining these non-aggregated variables, it becomes necessary to create the additional 
constraints that specify that there must be an aircraft available for each flight f £ C 
at some time. The reason we need these non-aggregated variables is to ensure that the 
necessary transfers of commodities occur at the flight level for continued flights. There 
are other methods of ensuring this transfer, but they involve adding a large number of 
side constraints to the formulation. Whereas, this method only requires the addition of 
\C\ constraints. 

The objective of the TFMRP is to minimize the total delay cost of flying all the 
required flights. Any flight may experience delay resulting from ground holding, decreas- 
ing speed while in the air, and selecting a route that is longer than the scheduled route. 
Moreover, a continued flight may also experience delay if there is no aircraft available for 
use at its departure time. 
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The objective function can be written in the following form: 

x d’xl.jt) + £ (» -<*/)«(«) +£«**&« + E 

{k,t,i=orig(k)} {( fc,t,/):/£C , {k,t,(i,j)&N(k)} 

k(f)=k > 

(4.1) 

The first term represents the cost of ground holding delay. The second term represents 
the cost of delay incurred by continued flights that were unable to depart on time as 
there were no available aircraft. The third term represents the cost of air delay due to 
speed reduction. Finally, the cost of delay caused by rerouting is obtained when the total 
cost of all scheduled air travel is subtracted from the total actual cost of air travel. The 
fourth term gives the total actual cost of all air travel and the fifth term is simply a 
constant representing the total cost of all scheduled air travel. 

The constraints are given below. 

(TFMRP) 

E - 12 ~ hi) = °» yie *. ( 4 - 2 ) 

12 W = 0 > Vk,t,i=dest(k),{ 4 . 3 ) 

= Dem k (t), Vk,t,i=dest{k),{4A) 
12 y/(0 + Xi Clic (t) - x$ Ctic (t- 1) - Xi A , ic {t~r{i))= 0, Vk,t,i=dest( k), (4.5) 

{/€C:fc=k'(/)} 

12 sfj*)" 12 Vf( t )- x % D ,i D ( t - 1 ) = Su Pi'( t )’ Vk,t,i=orig{k),{ 4.6) 

{f(i D ,})eN(k)} {/ec:k=k(/)} 

EE E < <*(*), Vi, t, (4.7) 

* 0:(t,j)e7V(fc)} 
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(4.8) 


£ y/(t) = 1, v/ s c, 

Oe?/} 

^•(<)>0, Vi, j\M, (4.9) 

y/(*)€{0,l}, V/, t. (4.10) 

Constraint (4.2) represents dynamic flow conservation for the sectors. There is a 
constraint for each sector node, i € <5, commodity fc, and time t. 

The next four constraints represent flow conservation at each of the airport nodes 
id, *b, *c, and ip. Constraint (4.3) forces flow conservation at node i A of airport i. At 
node %a, we sum over all the nodes that can arrive at airport i from some sector j of 
commodity k , 52{j-.(j,i A )€N(k)} x j,i A (* ~ t j,t A )- The time index is t - t hlA since this is the 

time that the flow leaves j if it is to arrive at i at time t. This flow must be equal to the 

flow out of node i A at time t. This flow goes to either node i B (where it will be removed 
from the network) or to node i c (where it will be transferred to another commodity). 

Constraint (4.4) forces flow conservation at node i B of airport i. The flow into node 
i B at time t equals the flow from node i A , x i AAB (t) plus any flow that is held on the 
ground from the previous time period, Xi gig (t - 1). This must equal the flow out which 
is Derrik(t ) plus x i BtiB (t). In other words, if flow arrives at node i B from node i A prior 
to the latest time that it is scheduled to arrive, it lands and then is held on the ground 
at no cost until the time that the flow can be removed with a nonzero value of Dem^t). 
Thus, every aircraft is removed from the network at its scheduled latest arrival time. 

Constraint (4.5) forces flow conservation at node ic- The flow into this node is of 
commodity k and the flow out of this node is a different commodity. The flow into node 
i c only comes from node i A , giving the term, E{h.i=destm x ti .<<?(* - r| d))- The term r ( i ) 
is the turnaround time at airport i. This is the time necessary in order to refuel and 
otherwise prepare the aircraft for the next flight. There is a delay arc at node i c which 
captures those aircraft that arrive at airport i and are cleaned and refueled before they 
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are needed to fly the next flight. There is no cost for using this delay arc, it simply 
represents an aircraft waiting at an airport for its next flight. Lastly, the flow out of 
node ic is given by £{/*=*'(/)} yf(t). This captures all the flights that are continued 
from commodity k flights that may depart at time t. 

Constraint (4.6) forces flow conservation at node id- At this node, all the flow leaving 
i D to some node j minus the flow into ip must equal the supply at airport i at time t for 
commodity k such that i = orig(k), which is denoted by Sup k (t). The flow leaving node 
i D is simply given by T,{j:(i D ,j)€N(k)} x i D , }(*■)■ This summation includes ground holding 
when j = ip- The flow into node ip is either from node ic or from ground holding in the 
previous time period. In the former case, the flow from ic passes through non-aggregated 
flight arcs. So this flow into node i o comes from all the flights of commodity k(f) — k 
that are continuations of previous flights. The ground holding value from the previous 
time period is given by x* D XD {t — 1). 

Constraint (4.7) captures the capacity restrictions. There is a capacity on the number 
of aircraft that can be within sector i at time t given by C,(f). To represent this in terms 
of the flow variables, we need to sum over all commodities and all arcs that represent 
travel in sector i of commodity k at time t. 

Constraint (4.8) forces the flights that are continued to occur. It states that for every 
flight that is a continuation of a previous flight, that flight must occur at some time, t, 
within a time window, 7). 

The remaining two constraint sets specify that the xj -(t) variables are nonnegative, 
and that the y/{t) variables are binary. Note that if y/(t) is binary, then x* ■(<) would be 
integer. 

The above formulation of the TFMRP is a multicommodity, integer variation of the 
minimum cost dynamic network flow problem. There are some important differences. 
First, the capacity constraint (4.7) is bundled over commodities, arcs and time peri- 
ods, not just over commodities. Second, the disaggregated variables yf(t), are not flow 
variables, and finally there are additional side constraints (4.8). 
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4.2 Modeling the Multiple Airline Problem 


The model proposed in the previous section can be used to solve the rerouting problem 
for a single airline. However, if we took the viewpoint of the FAA in which several airlines 
are occupying the airspace at the same time, then we need to modify the formulation 
slightly. The reason that the multiple airline problem is not the same concerns the 
continued flights. 

In the formulation of Section 4.1 we stipulate that a continued flight does not rely 
on a unique flight, rather the formulation guarantees that an aircraft of the correct 
commodity will be available before the continued flight can depart. Once we introduce 
multiple airlines, we need to make sure that not only is the incoming aircraft of the 
correct commodity, but also that it belongs to the correct airline. For instance, obviously 
an American Airlines aircraft could not continue a Delta Airlines flight. 

To do this we will need to redefine the commodities such that there is a unique com- 
modity distinguished by each origin-destination pair and by each airline. By redefining 
the commodities in this manner, we can now ensure that a continued flight will rely on 
an aircraft whose commodity corresponds to the correct airline. This would increase the 
number of commodities by a multiple equal to the number of airlines. 

We further need to consider the issue of fairness. It may be globally optimal to assign 
all the delay to a single airline, but of course, this solution is not acceptable. Thus, we 
would need to ensure that the delay is allocated in a fair manner across the airlines. This 
could be accomplished though modifying the Packing Formulation discussed in Section 

4.3.3 by adding constraints that guarantee that each airline receive no more than a given 
percentage of the total delay. 

4.3 The Lagrangian Generation Algorithm 

In this section we use Lagrangian relaxation of the formulation outlined in Section 4.1, 
randomized rounding, and a packing formulation to propose near-optimal solutions for 
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the TFMRP. The overall algorithm is outlined below. We then explain each step in 
detail. 

The Lagrangian Generation Algorithm 

1. (Lagrangian relaxation of the LP) Starting with the formulation of TFMRP, 
i.e., the problem of minimizing (4.1) subject to the constraints (4.2)— (4. 10) , we relax 
the capacity constraints (4.7) into the objective function with multipliers, A. We 
further relax the integrality constraints (4.10), i.e., we solve the relaxed problem as 
a linear program. We initialize the lower bound by solving the linear programming 
relaxation of the formulation for the TFMRP. The initial upper bound is infinity. 

2. (Solution of the relaxed problem) We solve the relaxed problem and obtain a 
potentially fractional solution y/(t ), 

3. (Randomized Rounding) We randomly round the variables y/{t) to zero-one 
solutions, and randomly decompose the flow into routes for flights. These routes 
are then added to a list of paths. 

4. (Packing formulation) We formulate and solve an integer packing problem, in 
which we are attempting to pack the elements of the list of paths into the capaci- 
tated airspace system. If a new solution is found, we update the upper bound. 

5. (Stopping criterion) If the upper and lower boimd are within a desired accuracy 
e, we stop. 

6. (Update of multipliers) We update the multipliers A and go to Step 2. 

In the following subsections we describe each of the steps of the algorithm in detail. 
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4.3.1 Lagrangian Techniques 

During Step 2 of the algorithm, we solve an uncapacitated multicommodity dynamic 
network flow problem as a linear program. Using the network flow solver of CPLEX, 
we solve the problem quickly (see Chapter 5). The main motivation of this step is to 
generate attractive routes for flights. 

We update the multipliers using the iterative approach of Everett (1963) as follows. 
We represent the capacity constraints (4.7) as Ax < b. Let a 3 be the jth row of the 
matrix A, and let bj be the right hand side value for this row. Let x k be the vector of 
solutions at the fcth iteration. Let A* be the value of the lagrange multiplier for the jth 
constraint at iteration k: 


If a'x* > bj then A} +1 = (1 + 6 k ) A*. 

If a < bj then A* +1 = (1 - 6 k )X k , 

( where the parameters 6 k are updated by the rule: 

If ( a' jX k — b } )(a' j x k ~ 1 — bj) > 0, then 6 k+l = ei6 k . 

If ( a!j X k — bj)(a' J x k ~ l — bj) < 0, then 6* +1 = e 2 S k . 

If ( a' 3 x k — bj){a'jX k ~ 1 - bj) = 0, then 8 k+l = 6 k . 

The values of and e 2 are fixed parameters where t\ > 1 and e 2 < 1. 

The motivation for this method is as follows. If a'x > b 3 then the solution x uses too 
much of the available amount of the jth resource; thus, we increase the lagrange multiplier 
in order to penalize the violation more. In this method, the lagrange multiplier would be 
increased by a factor of (1 + 6 k ). Likewise, if a'jX < bj then the solution x uses a feasible 
amount of the jth resource; thus, we decrease the lagrange multiplier. It is decreased by 
the factor (1 — S k ) as shown above. The amount of increase or decrease at each iteration 
i 
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is determined by 6*, called the step size, which is controlled at each iteration. 

The values of 6* are updated in the following manner. If a constraint is not satisfied 
iteration after iteration, then the step size is gradually increased based on the assumption 
that the value of A* may still be quite far from its optimal value. If the constraint 
fluctuates between feasibility and infeasibility, then the step size is reduced substantially 
based on the assumption that A(( has come close to its optimal value. It is interesting 
to note that updating the lagrange multipliers depends only upon whether or not the 
constraint was satisfied, not on the magnitude of the difference. 

Everett’s Method does not guarantee convergence. However, it has been shown to 
perform very well computationally; see Pugh (1993). 

4.3.2 Randomized Rounding Heuristic 

The objective of this step is to generate a rich set of paths for individual flights from 
the aggregated flow solutions. The motivation for using randomization is to generate 
a broader set of solutions. After completing the Lagrangian relaxation step of the La- 
grangian Generation Algorithm, we have a potentially fractional solution y/(t), 

Basically, this heuristic randomly walks through the network for every flight looking 
for a positive flow path. Starting at the departure airport of flight / and at the departure 
time of flight /, the heuristic randomly picks the next arc that has a positive flow on it. 
If this turns out to be a self-loop, then it remains at the airport for another time period 
and then makes another random decision about where to move in the next time period. 
Perhaps the next step places it at a new sector after it has completed the travel time for 
that arc. It then, once again, will randomly pick another arc that has some positive flow 
on it. It will walk through the network in this manner until it reaches its destination. 
With this in mind, we will now describe the process in full detail. 

We will create a list of paths as follows: 

P = {Pi,...,Pjv}, 
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where p t = {(s t (0), L(0)), . .., (^(n,), f,(n,))}, 


where s,(m) is the mth element of path p t , and t t (m) is the time that the flight arrives 
at s,(m), m = 0 , . . . ,n,. We define N as the total number of paths in the list and n, + 1 
as the total number of elements in path p t . 

We select the first element in each path Si(0) in a deterministic manner. For every 
noncontinued flight, we can create a path in which Si(0) is equal to node i D of the 
departure airport and t,(0) is equal to the scheduled departure time. For every continued 
flight, we select the earliest time such that y f (t) is non-zero and set this equal to one. 
Now we can create a path for each continued flight where $*(0) is equal to node i D of the 
departure airport and f,(0) is equal to the time at which y f (t) is equal to one. 

To build the rest of each path, we step through the network for every flight beginning 
at the node given by Sj(0) at the time t,(0). We will refer to the commodity of our flight 
as commodity k. Next, we randomly select from all arcs emanating from 5,(0) that have 
a positive flow value, i.e. z* ( o),j(^(0)) > 0. This may include the possibility of selecting 
the arc that represents ground delay. Let the arc that we randomly select by denoted by 
(si(0), j). We set 

s,(l) = j and U( 1) = *<(()) + f Si(0)J . 

We then need to decrease the flow value on the variable, x* -.(^(0)), by one. 

We continue in this manner until we reach the node representing node i J 4 of the 
destination airport. Since these flows respect the flow conservation constraints in the 
Lagrangian relaxation, there will always be flow out of a node that the heuristic reaches. 

There are a number of ways that we could set the probabilities used to select paths. 
Currently, we simply assign an equal probability to each node that has a positive flow. 
The rationale for this is simply to place a higher probability on obtaining alternative 
paths. Another possible method of randomizing would be to assign each arc a probability 
based on the flow on the arc. In particular, we could assign a probability P } to arc (a, j) 
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defined by 


T,{j:(a,j)eN(k)} x a,j{ t ) 

We could then use these probabilities to determine which arc to select at each step. 

After this heuristic is completed we have a set of path and time specifications that 
are added to a list of paths and used as input data for the Integer Packing Problem. 
Note that the heuristic only produces paths that satisfy flow conservation constraints. 
In particular, the capacity constraints (4.7), and as well as the airport constraints that 
handle continued flights (4.5), may not be satisfied. The path and time specifications 
will satisfy constraints (4.2), (4.3), (4.4), (4.6), (4.8), (4.9), and (4.10). By not forcing 
these paths to satisfy all the constraints in ( TFMRP ), we create, after a few iterations, 
a list of paths that has more flexibility. The Integer Packing Formulation will then select 
from this list a combination of paths that will satisfy all the constraints in (TFMRP). 

4.3.3 The Integer Packing Formulation 

The goal of the packing problem is to pack the paths generated by the randomized 
rounding heuristic into the air traffic system, so that all the necessary flights occur at 
correct (though not necessarily on time) times, and so that the capacity constraints are 
satisfied. 

A path, pi = (sj(0),t<(0», . . . , (Si(nj),t,(n i )), specifies the elements and the times of 
a given route. The elements are both airports and sectors. Obviously the first element, 
Sl ( 0) is the departure airport and the last element ^(n,) is the arrival airport. However, 
ground holding is represented in these paths. So, if path p, includes g units of ground 
holding then the elements s r (m ), m = 0, ...,g - 1, all represent the departure airport and 
the times U(m ), m = 0, ..., g - 1, are each separated by one time unit. Thus, time f,(0) is 
not the actual departure time, rather it is the time at which an aircraft becomes available 
to perform the flight represented by path p t . If this is a flight that is not continued, then 
t t (0) will be the same as the scheduled departure time. If the path represents a flight 
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that is continued, then t t ( 0) may be later than the scheduled departure time, since it is 
possible that there was no aircraft available for a continued flight to use at its departure 
time. 

The decision variables in the packing formulation are as follows: 



1, 

0, 


if path/time pair p, is used to fly flight / € T, 
otherwise. 


Let Z be the set of feasible combinations of path/time pairs and flights. 


Z- {(/,i): /€ C,p,(0) = orig(k(f)), p^n,) = dest(k(f)), U(0) > d f } U 
: f e f\C,Pi(0) - orig(k(f)), p^n,) = dest(k(f)), t t ( 0) = d f }. 

The objective of the packing problem is to minimize the cost of delay in the air as 
well as the ground holding cost of departing after the scheduled departure time. Let g 
once again be the number of units ground holding associated with path p t . Then the 
delay cost of path p x is given by 


c, = c a [U(ni) - ti(g)} + c?[ti(g-l) -**(())], 


which includes both the amount of time spent ground holding and flying. 
The objective function is as follows. 


N 


£ 


Ci 


Y, Z f,i + C 9 

{/■(/, t)€Z} 


N 


2 /.* X] df\ 

‘=1 {/:(/. t)€Z} /er J 


c a H. 


The first term captures the cost resulting from ground holding delay and from air 
travel. When combined with the final term, which is a fixed cost of scheduled air travel, 
we capture the cost of delay resulting from ground holding, decreasing speed while in the 
air, and selecting a route that is longer than the scheduled route. The only remaining 
delay occurs if there is no aircraft available for a continued flight to use at its departure 
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time. To capture this delay, we sum over all the times at which aircraft become available 
to perform the flight represented by path p t and subtract the sum of all the scheduled 
departure times. 

The constraint set is given by the following equations. 


( PP ) ]T z Jti < c 3 (t), Vj<=S,t, (4.11) 

{(f y i)€Z:3m\j=Si(Tn), 

£ Zf,i = 1, V/ € T, (4.12) 

{*:(/, «)€Z} 

y Zf t i< ^2 z p,i' i v/:, t, (4.13) 

{(/,<)ez:/ec, ')ez:fc(/')=fc, 

k'(/)=Mi (0)<t} (»V )+r(s,/ (n,i ))<t} 

2/,i€ {0,1}, V* = 1,...,JV. 

Constraints (4.11) represent the capacity constraints. They stipulate that for every 
sector ] and every time / . the sum over all the flights that are in this sector at time t 
must be less than or equal to the sector capacity. A flight is within a sector at time t if 
it entered the sector before time t, U(m) < t, and has not yet entered the next sector in 
its path before time t, t < ti(m + 1). 

Constraints (4.12) ensure that each flight will be assigned to exactly one route. 

Constraints (4.13) guarantee that a continued flight will not depart before a suitable 
aircraft has arrived for it to use. The left hand side of constraints (4.13) represent the 
number of continued flights whose preceding flight is of commodity k,k'(f) = k, and 
whose possible departure time is less than or equal to t. This number must be less than 
the number of flights of commodity k, k(f) = fc, that arrive before time t minus the 
turnaround time, r(sj/(rv)). 
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The final constraint set simply forces these variables to be binary. 


4.4 TFMP and MAGHP as Special Cases 

The air traffic flow management problem (TFMP) and the multi-airport ground- holding 
problem (MAGHP) can be seen as special cases of the rerouting problem. For instance, 
if we consider the TFMRP in which the choice of route through the network of sectors 
is fixed, then we have the TFMP. If we consider the TFMP in which the route consists 
of only a departure airport and an arrival airport, then we have the MAGHP. As both 
of these problems are of importance, it is reasonable to ask whether the Lagrangian 
Generation Algorithm could be used to efficiently solve these problems. 

4.4.1 The TFMP Dynamic Network Flow Formulation 

We shall start by describing how to model the TFMP as a multicommodity dynamic 
network. For TFMRP, we had variables of the following form. 


Jr / \ 

= number of flights of commodity k that depart from 
node i at time t and arrive at node j at time t + tij. 

The TFMP does not have as many choices for node j. In fact, for the TFMP, when 
the flow leaves node i there will be two options. It could be delayed at node i for one 
more time unit or it could proceed along its predefined path towards the next element 
in its path. If node i is one of the four airport nodes, then we have the exact same 
situation as we had for the rerouting case. If we define a path for each commodity as 
Rk = rfc,o> . . . , r kjlkl r klA }, where n k is the number of sectors, then we can modify 

the variables in the following manner. 



x k 0 (t) = number of flights of commodity k that are held at node r fc)i 
from time t to time t + l,for i = 0 
x^(t) = number of flights of commodity k that leave node r k>i 

at time t and arrive at rfc i+1 at time t + t k , for i = 0, ...,n k . 

The remaining variables are identical to the definition used for the rerouting problem. 

{ 1, if an aircraft scheduled to perform flight / arrives at time t, 

0, otherwise. 

By using these variables in (TFMRP), we obtain the following objective function for the 
traffic flow management problem. 

MIN Y. <? x Lin( t ) + £ (t-d f )y f (t)+ Y ca 4o(t)- 

{k,t,i=orig(k)} {t,f€C} {k,*,0<i<n*} 

The constraints are given by: 

+ = °. V M,* = i. 

<,i (t - t k nk ) - x k AtiB (t) - x k A lc {t) = 0, Vfc.t.t = dest(k), 

- 1 ) - 4 mb (*) = Dem k (t), \Jk,t,i=dest{k), 

Y yf( t ) + x ic,iA t ) ~ X L ~ x Lic ( t ~ r ( i ))= °> V M,* = dest(k), 

{/*=*'(/)} 
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{/*•«/» 

ds 

o 

II 

■ CO 

**-> 

> 

£ £ + (<'))<<?#), 
} {t' :t-t* <t* <t} 

Vi, t, 

5Z 0/(0 = !. 

V/ € C, 

4(0 > 0, 

Vi, j, k, t, 

y/(«)€{0,i} > 

V/, t. 

Note that many of the constraints are exactly the same as in ( TFMRP ). However, 
the majority of the constraints in (T FMRP) are due to the flow conservation constraints 
at sectors, (4.2) and from the capacity constraints (4.7). Both of these constraints are 

simplified for the TFMP due to the specification of the routes. 

We could now apply the 


Lagrangian Generation Algorithm to this dynamic formulation by relaxing the capacity 
constraints. The Integer Packing Formulation would be exactly the same as in Section 
4.3.3. 

4.4.2 The MAGHP Dynamic Network Flow Formulation 

For the ground holding problem we disregard the sectors and only consider the airports. 
So for the formulation we simply need the four node depiction of the situation at an 
airport as well as the corresponding constraints. We obtain the following objective. 

MIN £ 0 *?„,,„« +£ («-<*/) y,(t). 

{k,t,i=orig(k)} {t,fsC} 

The constraint set is given by the following. 
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VAr, t, i = dest.(k), 




“ X ) “***,*(0 =Dem k {t), Vk,t,i = dest(k), 


{/:*=*'(/)} 

4>,; D (*) + i, ix W - E y/(0 - *U(*-1) =Su ^( i )> Vk,t,i = orig(k), 
</:*-«/)> 

E s c -o(<). Vi„, i, 

{fc:t£>=ori0(/c)} 

E X L, <£>(*) + — Ci A (t), Vi A , t, 

{k:iA=dest(k)} 


E y/(0 = !. 

{ter,} 


v/eC, 


a£,(t) > 0, Vfc, t, (i, j) E N(k), 


y/ (t)e{o,i}, v/, t. 

This formulation differs from the last in that there are no sector flow conservation 
constraints and the capacity constraints have been separated into only departure and 
arrival capacity, given above as C lD {t) and C tA (£), respectively. Once again, we can now 
apply the Lagrangian Generation Algorithm to this formulation by relaxing the capacity 
constraints. 
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Chapter 5 


Computational Results for the 
Lagrangian Generation Algorithm 


In this chapter, we report on the computational performance of the Lagrangian Gener- 
ation Algorithm. In Sections 5.1, 5.2, and 5.3, the Lagrangian Generation Algorithm is 
applied to solve three instances of the air traffic flow management rerouting problem. 
These instances model a bad weather front passing through a portion of southwestern 
United States. This region consists of four airports, located at Denver (DEN), Phoenix 
(PHX), Las Vegas (LAS), and Salt Lake City (SLC). There are 42 sectors that lie in the 
vicinity of these four airports as shown in Figure 5-1. This data, that was provided by 
the FAA, includes the travel times for the sectors as well as the necessary turnaround 
times. Each of the three instances has a different capacity scenario determined by the 
movement of a weather front through the region. 


74 




< Figure 5-1: Sector map of the US with southwest region shown. 

5.1 Computations for Weather Scenario I 

For this instance, flight schedules for 71 flights between the 4 airports shown in Figure 
5-1 are extracted from a dataset provided by the FA A that covers an 8 hour time frame 
with 5 minute time intervals. We simulated a weather front passing from the northeast 
corner of this region to the southwest corner. We set the sector capacities according to 
a uniform distribution. 

During normal weather conditions, the sector capacities were generated using a uni- 
form distribution with a mean determined by the size of the sector and a standard devi- 
ation of one. Figure 5-2 shows the sector capacities during normal weather conditions. 
At the cusp of the storm, the capacities were generated using a uniform distribution with 
a mean of zero. All of the resulting negative values were set to zero. As the weather 
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Figure 5-2: Weather Scenario I at 8am representing normal operating conditions before 
the weather front has hit this region. 

front gradually passes through, the mean slowly increases until it reached the level dur- 
ing normal weather conditions. Figure 5-3 shows the different weather scenarios that 
evolve over time for this first capacity scenario as the weather front passes through the 
region. The shaded areas show the cusp of weather front with the numbers being the 
corresponding sector capacities. 

To achieve a lower bound on the solution, we solve the LP relaxation of the multi- 
commodity dynamic network formulation, which consists of 24,509 constraints and 61,912 
variables. Finding a solution to the LP requires 181 seconds. The solution is highly non- 
integral with an objective value of 2498.5. 

We ran the Lagrangian Generation Algorithm setting the parameter values as ej = 2 
and 62 = 0.33. The starting values for A ® and 8 ® were set to 10 and 0.8 respectively. 
We did not perform extensive trials to determine the best starting values for A and 
however, a few different settings were tried. The results presented reflect the stated 
starting conditions, which converged in the shortest number of iterations. The Lagrangian 
relaxation solves a network problem with 15,279 nodes, 54,427 arcs and 26 of the side 
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Figure 5-3: Weather Scenario I at: (a) 9:45am (b) 10:35am (c) 11:25am (d) 12:15pm. 

constraints (4.8). The size of the integer packing problem grows at each iteration as 
the size of the list of paths increases. At the final iteration, the formulation consists of 
2209 constraints and 145 variables, which reduces to 424 constraints and 125 variables 
by using some presolving routines in CPLEX. Table 5.1 tracks the performance of the 
Lagrangian Generation Algorithm as it steps through the algorithm. All of times are 
given in seconds. 

The total amount of time needed to solve this problem, including the subproblem 
times for the Lagrangian relaxation and the integer packing problem given in Table 5.1, 
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Lagrangian 

Packing 

Iteration 

Time 

Obj Value 

No. Infeas. 

Presolve Time 

IP Time 

Obj Value 

0 

25.01 

-113062 

89 

0.03 


Inf. 

1 

24.22 

-12572 

65 

0.03 


Inf. 

2 

24.70 

1321.46 

77 

0.07 

- 

Inf. 

3 

24.42 

2073.94 

54 

0.12 

- 

Inf. 

4 

24.58 

2248.13 

64 

0.18 

- 

Inf. 

5 

24.40 

2324.27 

68 

0.18 

- 

Inf. 

6 

24.65 

2286.42 

61 

0.18 

- 

Inf. 

7 

24.67 

2377.25 

66 

0.18 

15.26 

2509 


Table 5.1: Computational Results for weather scenario I. 


was 330 seconds. The solution value found, 2509, is within 0.4% of the lower bound. 
The total amount of ground holding imposed upon this solution was 810 minutes and 
the total amount of airholding was 15 minutes. 

Numerous routes were used for each particular commodity. As an example, we will 
consider one of these commodities, Las Vegas to Phoenix, in detail. Three different routes 
were used over the course of the day to fly flights of this commodity. Two routes are 
1 shown and labeled in Figure 5-4. Most of the time when this commodity is flown, route 



Figure 5-4: Two routes used at different times to fly from Las Vegas to Phoenix. 

1 is used and ground delay alone is assigned to avoid any anticipated capacity problems 

1 


78 





while in the air. However, we will examine two flights that use route 2. The important 
times for these two flights are depicted in Figure 5-5 and explained below. 


Flight l 


45 mins 


actual 

departure 

i2:35 


I hr, 35 mins 
air travel 


actual 

arrival 

2:10 


ground holding* 

1 


r t 

11:50 

1 

12:50 

j 

scheduled 

scheduled 


departure 

arrival 



Flight 2 

actual 

arrival of actual . actual 

aircraft departure 1 ‘ arrival 
11:25 14 :35 airtravel 1^20 


11:00 H ;20 12:20 

scheduled scheduled scheduled 

arrival of departure arrival 

aircraft 


Figure 5-5: Time lines of two flights that use route 2 on Figure 5-4. 

At 11:50am a flight is scheduled to depart from Las Vegas. This flight is scheduled 
to arrive at Phoenix one hour later. The flight is held on the ground for 45 minutes and 
actually departs at 12:35pm. It reaches Phoenix at 2:10pm making the total amount 
of time spent traveling equal to 1 hour and 35 minutes, which is 35 minutes more than 
scheduled. The total delay of 80 minutes resulted from 45 minutes of ground holding and 
35 minutes of rerouting delay. 

At 11:20am, a continued flight is scheduled to depart from Las Vegas. It is relying on 
an aircraft that is scheduled to arrive at Las Vegas at 11:00am. However, this incoming 
flight experiences delay and does not arrive at Las Vegas until 11:25am. The aircraft 
is immediately refueled and boarded, and the continued flight actually departs at time 
11:35am, 15 minutes after the scheduled departure time. There is no ground holding 
assigned to this flight. The flight travels along route 2 and is further delayed in terms 
of airspeed reduction. When the flight reaches the shaded sector in Figure 5-4, its speed 
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4 


is reduced such that the travel time through this sector is increased by 10 minutes. It 
reaches Phoenix at 1:20pm, a total of 1 hour after the scheduled arrival time. This 
total delay of 60 minutes resulted from 10 minutes of airspeed reduction, 35 minutes of 
rerouting delay and 15 minutes of departure delay due to the late incoming aircraft. 

5.2 Computations for Weather Scenario II 

For this instance, the same flight schedules for the 71 flights between the 4 airports 
shown in Figure 5-1 were used. However, to simulate the weather front passing from the 
northeast corner of this region to the southwest corner, the sector capacities were set 
deterministically. 

During normal weather conditions, we fixed the capacities according to the size of the 
sector. At the cusp of the storm, the sector capacities were set to zero. One hour later, 
as the storm front moves along, those sectors that had zero capacity during the last hour, 
now have a slightly increased capacity of one. The sector capacities would continue to 
increase hourly until they have resinned the level of normal weather conditions. Figure 
5-6 shows the different sector capacities that evolve over time for the second capacity 
scenario as the weather front passes through the region. Figure 5-6 (a) through Figure 5- 
6 (f) show the sector capacities between the times 8:20 to 9:20am, 9:25 to 10:25am, 10:30 
to 11:30am, 11:35am to 12:35pm, 12:40 to 1:40pm, and 1:45 to 2:45pm, respectively. The 
shaded areas show the cusp of the weather front in which the sector capacities are zero. 
The front appears to move along gradually spending one hour before each progression. 

To achieve a lower bound on the solution, we solve the LP relaxation of the multicom- 
modity dynamic network formulation. The size of the formulation is not affected by the 
change in the weather capacity scenario. Thus, the number of constraints and variables 
is the same as specified for weather scenario I. Solving the LP requires 59 seconds and 
gives a solution that is highly non-integral with an objective value of 2387. 


< 
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Figure 5-6: Weather Scenario II: (a) 8:20 to 9:20 am, (b) 9:25 to 10:25am, (c) 10:30 to 
11:30am, (d) 11:35am to 12:35pm, (e) 12:40 to 1:40pm, (f) 1:45 to 2:45pm. 
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Lagrangian 

Packing 

Iteration 

Time 

Obj Value 

No. Infeas. 

Presolve Time 

IP Time 

Obj Value 

0 

24.72 

-107072 

72 

0.12 

- 

Inf. 

1 

24.63 

-11148 

60 

0.10 

- 

Inf. 

2 

24.72 

1618.38 

20 

0.42 

17.51 

2418 


Table 5.2: Computational Results for weather scenario II. 


We ran the Lagrangian Generation Algorithm setting the parameter values as q = 2 
and e 2 = 0.33. The starting values for and 8 ® were set to 10 and 0.8 respectively. 
Table 5.2 tracks the performance of the Lagrangian Generation Algorithm as it steps 
through the algorithm. The problem sizes are the same as for weather scenario I. All of 
times are given in seconds. 

The total amount of time needed to solve this problem, including the subproblem 
times for the Lagrangian relaxation and the integer packing problem given in Table 5.2, 
was 116 seconds. The solution value found, 2418, is within 1.2% of the lower bound. 
The total amount of ground holding imposed upon this solution was 480 minutes and no 
airholding was used. 



Lagrangian 

Packing 

Iteration 

Time 

Obj Value 

No. Infeas. 

Presolve Time 

IP Time 

msraflBEi 

0 

24.72 

-107072 

72 

0.12 

- 

Inf. 

1 

24.63 

-11148 

60 

0.10 

- 

Inf. 

2 

24.72 

1618.38 

20 

0.42 

17.51 

2418 

3 

26.18 

2104.20 

25 

0.48 

17.42 

2408 

4 

25.98 

2231.75 

32 

0.53 

18.02 

2396 

5 

25.98 

2295.00 

18 

0.53 

19.10 

2393 

6 

25.98 

2317.74 

19 

0.57 

19.35 

2392 

7 

25.43 

2345.69 

21 

0.62 

18.46 

2389 j 

8 

25.88 

2327.27 

21 

0.65 

19.56 

2389 

9 

25.82 

2360.46 

22 

0.65 

19.89 

2389 

10 

25.88 

2304.59 

22 

0.68 

19.99 

2389 

11 

25.78 

2367.44 

10 

0.72 

20.03 

2389 

12 

26.55 

2375.27 

12 

0.68 

20.36 

2389 

13 

25.93 

2378.53 

14 

0.72 

21.58 

2389 

14 

26.45 

2379.18 

21 

0.67 

21.60 

2389 


Table 5.3: Computational Results for weather scenario II for 15 iterations. 
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We ran the algorithm for one hundred iterations, without the e stopping condition, 
to see if we could find a solution which is even better than the one found above. The 
results from the first fifteen iterations are given in the Table 5.3. The remaining eight- 
five iterations did not find a solution with a lower solution value less than the value at 
the fifteenth iteration, and were thus not including in the table. Within 7 iterations we 
generate the best solution that we were able to find. This value is 2389, which is within 
0.08 % of the lower boimd of 2387. 



Figure 5-7: Three routes used to fly between Salt Lake City and Phoenix. 

Numerous routes were used for each particular commodity. In Figure 5-7 we look at 
a few of the routes used to fly between Phoenix and Salt Lake City. Route 1 travels from 
Salt Lake City heading towards Phoenix. This flight is scheduled to depart at 9:30am, 
but is held on the groimd until 11:30am, incurring a two hour ground delay. At that 
point, the flight departs from Salt Lake City and follows a reasonably direct route, route 
1, to Phoenix, basically trailing the storm front. 

Routes 2 and 3 go north from Phoenix to Salt Lake City. The flight that travels along 
route 2 is scheduled to depart at 9:25am, but is held on the ground for 1 hour. It then 
departs and follows a circuitous route due to the limited sector capacities. The flight 
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that travels along route 3 is scheduled to depart at 11:30am and only suffers a 20 minute 
ground delay. This flight passes through its route immediately after the weather front. 
The capacities are still quite limited, however, forcing this flight to significantly deviate 
from the shortest path. In order to appreciate how this affects travel time, we note that 
route 1 requires 1 hour and 55 minutes of flying time, route 2 requires 2 hours and 25 
minutes of flying time and route 3 requires 2 hours and 35 minutes of flying time. 

5.3 Computations for Weather Scenario III 

The final weather scenario that we tested uses the same geographic area as shown in 
Figure 5-1. However, we increased the number of flights to 200 flights and scaled the 
normal weather sector capacities accordingly. To simulate the weather front passing from 
the northeast corner of this region to the southwest corner, we set the sector capacities 
deterministically. 

During normal weather conditions, we fixed the sector capacities according to the 
size of the sector. At the cusp of the storm, the sector capacities were set to zero. The 
front once again gradually moves along spending forty minutes, before each progression. 
As the cusp of the storm front moves through the region, the available sector capacity 
increases by two each forty minutes. So this weather front moves more quickly and does 
not leave such bad conditions behind it, as does the previous weather scenarios. Had we 
kept the capacity levels at that of either the two previous scenarios, then the problem 
would have been infeasible, meaning that there would have been no way to complete all 
of the 200 flights during the time frame without cancelling some flights. Figure 5-8 shows 
the different weather scenarios that evolve over time for the third capacity scenario as 
the weather front passes through the region. Figure 5-8 (a) through Figure 5-8 (f) show 
the sector capacities between the times 8:20 to 9:00am, 9:05 to 9:45am, 9:50 to 10:30am, 
10:35 to 11:15am, 11:20am to 12:00pm, and 12:05 to 12:45pm, respectively. The shaded 
areas show the cusp of the weather front in which the sector capacities are set to zero. 
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Figure 5-8: Weather Scenario III for 200 flights: (a) 8:20 to 9:00am, (b) 9:05 to 9:45am, 
(c) 9:50 to 10:30am, (d) 10:35 to 11:15am, (e) 11:20am to 12:00pm, (f) 12:05 to 12:45pm. 
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To achieve a lower bound on the solution, we solve the LP relaxation of the multi- 
commodity dynamic network formulation, which consists of 25,881 constraints and 66,489 
variables. Solving the LP requires 86 seconds and gives a solution that is highly non- 
integral with an objective value of 6513.5. 

We again ran the Lagrangian Generation Algorithm with the same starting values 
as before. Table 5.4 tracks the performance of the Lagrangian Generation Algorithm 
as it steps through the algorithm. All of the times given in Table 5.4 are in seconds. 
The Lagrangian relaxation solves a network problem with 16,219 nodes, 57,294 arcs and 
138 of the side constraints (4.8). The size of the integer packing problem grows at each 
iteration as the size of the list of paths increases. At the final iteration, the formulation 
consists of 2394 constraints and 1197 variables, which reduces to 879 constraints and 
1078 variables by using some presolving routines in CPLEX. 



Lagrangian 

Packing 

Iteration 

Time 

Obj Value 

No. Infeas. 

Presolve Time 

IP Time 

Obj Value 

0 

36.52 

-283953 

91 

0.19 

- 

Inf. 

1 

36.49 

-26652 

54 

0.26 

- 

Inf. 

2 

36.31 

5772.76 

39 

0.81 

- 

Inf. 

3 

36.56 

6152.06 

44 

0.83 

18.11 

6574 


Table 5.4: Computational Results for 200 flight dataset. 


The total amount of time needed to solve this problem, including the subproblem 
times for the Lagrangian relaxation and the integer packing problem given in Table 5.4, 
was 169 seconds. The solution value found, 6574, is within 0.92% of the lower bound. 
The total amount of ground holding imposed upon this solution was 700 minutes and no 
air holding was used. 

Once again, we ram this problem for more iterations to see if we could obtain a better 
solution, even though we already found a solution that is well within the tolerance. Table 
5.5 gives results that show all the solutions that we generated. Of these, the best solution, 
6520, is within 0.09% of the lower bound, 6513.5. 

Again, we can see how well the Lagrangian Generation Algorithm performs in terms 
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Lagrangian 

Packing 

Iteration 

Time 

Obj Value 

No. Infeas. 

Presolve Time 

IP Time 

Obj Value 

0 

36.10 

-283953 

91 

0.19 

- 

Inf. 

1 

36.16 

-26652 

54 

0.26 

- 

Inf. 

2 

36.14 

5772.76 

39 

0.81 

- 

Inf. 

3 

36.26 

6152.06 

44 

0.83 

18.11 

6574 

4 

36.08 

6291.80 

42 

0.92 

18.78 

6545 

5 

36.03 

6327.60 

39 

0.94 

18.63 

6529 

6 

35.97 

6307.54 

42 

1.00 

19.00 

6527 

7 

35.83 

6370.08 

37 

1.04 

19.12 

6523 

8 

36.51 

6453.82 

36 

1.08 

19.58 

6520 

9 

36.26 

6458.95 

47 

1.09 

19.64 

6520 

10 

35.90 

6319.08 

41 

1.10 

20.01 

6520 

11 

36.22 

6380.18 

36 

1.10 

20.82 

6520 


Table 5.5: Computational Results for 200 flight dataset for 12 iterations. 


of speed and accuracy of solutions. This last dataset is quite large, but the size of the 
multicommodity dynamic network flow formulation is not drastically increased by an in- 
crease in the number of flights. Only the number of non-aggregated variables is increased; 
the number of aggregated variables remains the same. The only set of constraints that is 
affected by this increase is constraint (4.8), which comprises a very small amount of the 
total number of constraints. 
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Chapter 6 


Conclusions 


In this thesis, we have presented what we believe to be interesting and practical ap- 
proaches for solving the Air Traffic Flow Management Problem and the Air Traffic Flow 
Management Rerouting Problem. In Section 6.1, we summarize the contributions of the 
research presented in this thesis. In Section 6.2, we discuss some remaining issues of 
concern to the FA A. 


6.1 Contributions 

The Air Traffic Flow Management Problem determines how to allocate ground holds to 
flights as well as how to control the enroute speed of an aircraft while taking into account 
all the capacitated elements in the system (arrival, departure and sector capacities). We 
built an integer programming formulation to minimize delay costs. The formulation is 
quite strong as some of the proposed inequalities are facet defining for the convex hull 
of solutions. We addressed the complexity of the TFMP and showed that it is NP-hard. 
The formulation easily extends to incorporate the dependence of airport runway capacity 
of departures and arrivals, hub connectivity, and banking of flights. When specialized for 
the multiple airport ground holding problem, we proved that the LP relaxation bound of 
our formulation is at least as strong as others proposed in the literature. The solutions 


! 


88 



of the LP relaxation of the TFMP were almost always integral, so there was no need 
to branch and bound. In essence, our formulations reduced the problem to efficiently 
solving large scale linear programming problems. As a result, the computation times were 
reasonably small for large scale, realistic size problems involving thousands of flights. 
Short computational times and integrality properties are particularly important, since 
these models are intended to be used on-line and solved repeatedly during a day. 

The Air Traffic Flow Management Rerouting Problem (TFMRP) determines how to 
reroute flights through different flight paths in order to reach their destinations if the 
current routes pass through a region that is unusable as a result of poor weather con- 
ditions. This is the first research that has taken a global look at rerouting. We not 
only determined the best routes for the aircraft to follow, but also the amount of ground 
holding, and the amount of speed adjustment while taking into consideration that the 
entire national airspace system is operating under a strict capacity restriction. To address 
this problem, we presented an integrated mathematical programming approach, which 
utilizes several methodologies, for the problem of minimizing delay costs. In order to 
address the high dimensionality, we began by presenting an aggregate model, in which 
the problem is formulated as a dynamic, multicommodity, integer network flow prob- 
lem with certain side constraints. Using Lagrangian relaxation, we generated aggregate 
flows that are decomposed into a collection of flight paths for individual aircraft using a 
randomized rounding heuristic. This collection of paths was then used in a packing inte- 
ger programming formulation, the solution of which generates feasible and near-optimal 
routes for individual flights. The overall algorithm, termed the Lagrangian Generation 
Algorithm, was used to solve real problems in the southwestern portion of United States. 
The solutions returned by our algorithm were within 1% of the corresponding lower 
bounds. However, there are several remaining issues that need to be addressed before 
this approach can be implemented. In Section 6.2, we discuss some of these issues. 

In the course of this research, we obtained some general insights that may have 
wider applicability. First, we have shown through our formulation of the TFMP, that 
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by redefining traditional time assignment variables as “by time t” rather than “at time 
t " , some of the proposed inequalities are likely to be facet defining for the convex hull 
of solutions, which in turn leads to considerable computational improvement in terms of 
achieving integrality. We believe that this redefinition may potentially be beneficial for 
other scheduling problems where we are assigning items to specific time intervals. 

Second, the algorithmic design of the Lagrangian Generation Algorithm could be used 
in other problem contexts. The idea of extracting solutions using randomization and 
combining these solutions using integer programming may be useful in other problems 
as well. 

Although we have presented our formulations in the context of air traffic control, we 
envision other applications of our models in any area in which goods are dynamically flow- 
ing through a system with several types of capacitated elements such as manufactming, 
telecommunications and ground transportation systems. 

6.2 Other FA A Issues 

For several years, the FAA has been operating an Air Traffic Control System Com- 
mand Center (ATCSCC) in Washington, D.C. This center is equipped with outstanding 
information-gathering capabilities that dynamically keeps track of all the information 
about capacities, flight information, weather, etc. As we have mentioned earlier, the FAA 
uses a computerized procedure to allocate ground holding delays based on a first-come- 
first-serve rule. We believe that the optimization-based approaches that are described in 
this thesis are well suited to be the optimization “brain” for this system. However, there 
are important issues that need to be addressed before applying an optimization based 
approach in a real world environment. 

1 Interaction with airlines. 

After the ground delays are issued, the airlines have the opportunity to propose 
modifications to these delays through a cancellation and substitution process. It 
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would be important to elucidate the effects of this interaction. 

2 Passenger and flight crew connections. 

There are more connectivity requirements than just the aircraft connectivity. In 
particular, there are the connections of the passengers and the flight crew. On 
occasion delay will be incurred if many of the passengers have not yet arrived due 
to delay encountered during the first leg of their journey. Moreover, if the delay 
involves the flight crew, either an alternative flight crew must be reassigned or the 
flight will experience a delay in waiting for the crew to arrive. 

3 Mechanical Difficulties. 

Flights are often further delayed by mechanical failures of the aircraft. Research 
could be done that allows for the possibility of randomly occurring delay from 
equipment failure. 

4 Flight cancellations. 

When weather conditions become bad enough to warrant extreme measures, flights 
are often cancelled or forced to land at alternative destinations. This situation 
causes many difficulties throughout the air traffic system. However, it is a viable 
option that has not yet been considered. 

5 Dynamic updating of decisions. 

Both ground and enroute delays are determined simultaneously several hours before 
a flight leaves. In practice, however, enroute delays are not known until after the 
aircraft is in the air. Clearly more research is needed to clarify the implications of 
considering enroute delays on a much shorter time scale and to demonstrate how 
to effectively update the previous solution to incorporate any new information. 

6 Stochastic modeling. 

The model presented in this paper assumes a deterministic environment. More re- 
search is needed to account for stochasticities inherent in a system that is strongly 
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dependent upon weather conditions. 
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Appendix A 


On the Polyhedral Relationships 
Between Ground Holding 
Formulations 


We intend to establish Proposition 1 . Since I P'maghpi IPvbo I Ptp af e valid integer 
programming formulations, it is clear that IP'maghp = IPvbo = IPtp ■ Moreover, since 
the IP is more restrictive than its relaxation, IP'maghp — p maghp- 

To show the relationship P'maghp Q P tp we wil1 start with a feasible point in 
P'maghp > ^/t, and show that this is indeed feasible to Ptp- The first two constraints 
and the last constraint are identical in the two models. So what remains to be shown 
is that any point, ~Zft, that satisfies the third constraint of P'maghp a ls ° satisfy the 
third constraint of Ptp- So 


Y (zft - - Y ( z ft' - 1 ) 

t€T f a ,t<r t'€T?i,t'<T—s f i—(rf —df) 

= Zfr Z S' ,T—Sj> — (rj—dj) z f,t+(rf— df) Z f,t — Sji — 0 

So the point, z/ f , satisfies the third constraint of P T p a nd all of the constraints of 
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Ptp hold. Thus, the point z/ t does indeed lie in the polyhedron Ppp- This establishes 
the relationship P'maghp — Ptp- 

Now we need to prove the relationship P TP C Pvbo- To show this we will start with 
a feasible point in Ptp, Zft, and show that this is indeed feasible to Pvbo- Once again, 
the first two constraints and the last constraint are identical in the two models. So what 
remains to be shown is that any point, Zf t , that satisfies the third constraint of P TP will 
also satisfy the third and fourth constraints of Pvbo- So 

9f = £f(«/t-T/ lt _i)-r / 

t£Tj 

9f ~ df + Sf> + ry — gj 

— X] ~ 1) — r f — df + Sf> + ry - ^2 t i^ft - Zf,t- 1 ) - Tf 

terp terp 

= - ... - Zf >r > f +T f ,-1 + r f + Tf ~ r f -d f + Sy + r f , 

+z ft r f + • ■ • + z fT{+ f f _ 1 -r f -Tf + r f 

= -Zf'y } - ... - ~Z f> t r'j+T fi — \ + T f ~ df + Sy + Ty + Z/, r/ + . . . + Z frf+ f f _ x ~ Tf 
— —( r f + Ty — 1 — (r f — Sy — ( Tf — df)) + 1) + Ty — df + Sy + Ty 

+(Tf + rf-l-r f + l)-T f 

= — (ry + Ty — Tf + Sy + rf — df) + Ty — df + Sy + Ty + 0 

= —rfi — Sy + df — df + Sy + Ty 
= 0 

where Tf is maximum amount of time that flight / may arrive late. So all of the con- 
straints hold and the point Zf t does indeed lie in the polyhedron Pvbo- This establishes 
the relationship Ptp C Pvbo- ■ 
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Appendix B 


On the Non-integrality of the 
Polyhedron PmaGHP 


In this section we prove Theorem la, i.e. , the polyhedron Pmaghp is n °t integral, by 
providing the following example which has a fractional extreme point. Consider the case 
in which there are two flights arriving and being continued by two flights departing from 
a given airport during a restricted time window. The data of the problem is as follows: 

|/C| = 1, T = {1,2, 3, 4}, C = {(1,1), (2,2)}, 

i.e., the arriving flight i is continued by departing flight i. The turnaround times are 
Sl = 0, s 2 = 1. The time windows are: 

7? = {1,2}, T 2 Q = {1,2}, Tf = {1,2}, T 2 d = { 2,3}. 

Notice that flight 2 can only depart during time slots 2 and 3 since the turnaround time 
for the second flight is 1. The decision variables are: 

2/11 , yi2, 2/22 , 2/23) 2 11> 2 12) z 21 , z 22> 
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with the interpretation that y,j = 1 if flight i departs by time j and z tJ = 1 if flight i 
arrives by time j. Because of the time windows, we know that 

2/13 = 1. 2/24 = 1 r. Z\3 ~ 1) 223 = 1- 


The capacities are: 

D{ 1) = D( 2) = D( 3) = 1, A(l) = >1(2) = >1(3) = 1. 

The resulting formulation ( MAGHP ) is: 

2/11 < 1 , 211 +221 < 1 , 2/12 

2/12 - 2/11 + 2/22 < 1, 212 - 211 + Z22 - 221 < 1, Zl2 

1 - 212 + 1 ~ 2 22 < 1, 2/11 

2/12 


(2/11 ) 2/12, 2/22, 2/23 1 2n, Z12, Z21 > 222 / 

( 1 , 1 , 0 , 1 , 1 , - 1 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 )' 

and 


1 - 2/12 + 2/23 - 2/22 < 1, 


Letting 

rc — 

6 = 


— 2/11 > 0, 2/23 - 2/22 > 0, 

- 211 > 0, Z22 - 221 > 0 , 

- 2n < 0 , 2/22 - 221 < 0 , 

— Z12 < 0 , 1/23 — 222 < 0 
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yn Vl2 2/22 2/23 Z U z \2 z 2 \ z 22 

/lOOOOOOO^ 
-1 1 1 0 0 0 0 0 

0 -1 -1 1 0 0 0 0 

0 0 0 0 1 0 1 0 

0 0 0 o-ll-ll 

0 0 0 0 0 -1 0 -1 

1 -1 0 0 0 0 0 0 

001-10000 
0 0 0 0 1 -1 0 0 

0 0 0 0 0 0 1 -1 

1 0 0 0 -1 0 0 0 

0 1 0 0 0 -1 0 0 

0 0 1 0 0 0 -1 0 

^ 0 0 0 1 0 0 0 -1 j 


the feasible space can be written as Ax < b. 

Notice that matrix A is not totally unimodular since the submatrix consisting of the 
columns corresponding to the variables yi2 ■ ?/22 ■ ~ 12 ■ and 221 and the third, fifth, twelfth 
and thirteenth rows: yi2 y22 2l2 * 21 

/ —1 —1 0 0 ^ 

001-1 
10-10 
v 0 1 0-1 J 

has determinant of 2 . The objective function 


Min 2 y n - 4 y ]2 + 2 y 2 2 - 6 1/23 - 3 z n + 6 212 - 3 2 2 i + 6 *22 
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gives an optimal solution of 


Vn = 0 2/22 = 0 z n = \ z 2 \ = \ 

OJ 1 „ _ 1 ^ 1 ^ _ 1 

2/l 2 2 2^23 2 *^12 2 ^22 2 

that shows that the polyhedron Pmaghp is not integral. Furthermore, this is the objective 
function that is obtained when we let cj = 1, c“ = 3 for all / € T. So, even with the 
restriction that c g = cj and c a = cj for all / € T, the polyhedron Pmaghp is not integral. 
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Appendix C 


Facet Defining Constraint Proofs 


In this section we analyze the polyhedral structure of the ccm,v(I P maghp) and provide 
the proof of the first half of Theorem lb that establishes which constraints are facets 
of conv(I Pmaghp)- The proof of the second half of Theorem lb concerning problem 
(TFMP) is similar, but more algebraically involved. We first show that the constraint 

E (y/t-yu-i)<Dk(t), vk e jc,t e f 

is not a facet of conv (I P maghp) by constructing a counterexample with two flights, one 
arriving at airport k and one departing from airport k, three time periods and D(t) = 1 , 
A(t) = 1 . Then only the variables y n ,y x2 , 2/13, Z11, 212, and z 13 are defined. The complete 
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set of feasible solutions to IPmaghp is given by: 


y n 

0 

3/12 

0 

!/13 

0 

211 

0 

212 

0 

213 
0 ^ 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1 

1 

0 

0 

0 

1 

1 

1 

0 

0 

1 

0 

0 

1 

0 

0 

1 

0 

1 

1 

0 

0 

1 

1 

1 

1 

0 

1 

1 

0 

1 

1 

0 

1 

1 

1 

1 

1 

V 1 

1 

1 

1 

1 

1 J 


In this case, dim(I Pmaghp) — 5 which can be determined by checking the rank of 
the matrix of solutions. We define the set 


H t = {( y , z) e IPmaghp ■ ~ Vu-i) = 1}> f° r some t € ?■ 

{f-terf} 

Then, 

H z = {(0,0, 1,0, 0,1), (0,0, 1,0, 1,1), (0,0, 1,1, 1,1)}. 

In this case, the maximum number of affinely independent points in i/3 is less than the 
dim(I Pmaghp) — 1 - We conclude that the constraint H{/ : t G r*}(y/t — — Dk{t),Vk € 
fC,t € T is not a facet. The same result can be checked in a similar manner for the 
constraint < A k (t), Vk,t. D 

For ease of exposition we consider instances of (M AG HP) such that 

• \T/\ is that same for all / and therefore D = maxf\Tf \ = \Tj\, 

• s f = o,v/ e P, 
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• A k (t),D k (t) > 1 


We consider an instance of {MAG HP) with |JF| flights in which \C\ (< 1^1) of these 
flights are continued. These flights were arranged such that the first \C\ flights are 
continued by flights \C\ + 1, . . . , 2\C\ < \T\, with flight 1 being followed by flight \C\ + 1, 
flight 2 being followed by flight \C\ + 2, and so on. 

We first determine dim(I Pmaghp ) by constructing the following matrices of solutions, 
in which each row represents a solution to {MAG HP), (see Figures C-l and C-2). The 
rows of these matrices are affinely independent and there are 2\T\D + 1 such rows. 
So, we have exhibited 2\P\D + 1 affinely independent points in IPmaghp and thus, 
dim(I Pmaghp) — 2\P\D. 

We next consider the set 

G ft = {(y,z) € IPmaghp ■ Vft - Vu-i = 0}, for some f eF,t<=T. 

If / € {1, . . . , |C|} then there are four distinct solutions from the matrices of Figures 
C-l and C-2 which do not belong to G ft . For each of these rows, replace the 0 in the 
yf t - 1 column with an 1. 

If / € {|C| + 1, . . . ,2|C|} then there are two distinct solutions from Figures C-l and 
C-2 which do not belong to G/ t . For each of these rows, replace the 1 in the y/. t column 
with a 0. 

If / e {2|C| + 1, . . . , |^|} then there are two unique solutions from Figures C-l and 
C-2 which do not belong to G/ ( . For each of these rows, replace the 0 in the y/,t_i column 
with a 1. 

For all of these cases, we have constructed a matrix with \T\ D affinely independent 
rows, proving that dim(Gft) > \P\ D — 1. Since Gft is a proper face of IPmaghp > we 
know that dim(Gft) < dim(I P maghp)- So, dim(Gft ) = — 1 and thus, Gft is a 
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facet of I Pm ag hp ■ □ 

We next consider the set 


Kft = {(j/< z) G IPmaghp '■ Zft ~ 2 /,t-i — 0}, for some / G T,t G T . 

If / G {1, . . . , |C | } then there are three distinct solutions from the matrices of Figures 
C-l and C-2 which do not belong to G j t . For each of these rows, replace the 1 in the yj t 
column with a 0. 

If/G {|C| + 1,...,|F|} then there is only one distinct solution from Figures C-l and 
C-2 which does not belong to G/ t , so remove this row. 

For each of these cases, we have constructed a matrix with \J-\D affinely independent 
rows, proving that dim(Kf t ) > \P\ D — 1. Since Kf t is a proper face of IPmaghp, we 
know that dim(Kf t ) < dim(I Pm ag h p) • So, dim(Kj t ) = \T\D — 1 and thus, Kft is a 
facet of IPmaghp- □ 

We next consider the set 


Mft = {(y,z) G / Pmaghp ■ Zft - y fj—(r f—df) = o}, for some f e T,t <eT. 

For all / G { 1 , . . . , \P\ } there are t — Tj + 1 distinct solutions from the matrices of Figures 
C-l and C-2 which do not belong to M fl . For each of these rows replace the 0’s in the 
columns corresponding to Zf t >,t < t' < T / with l’s. Tf and T/ are the last possible and 
the earliest possible times that flight / could arrive, respectively. 

The remaining matrix will have \IF\D affinely independent rows, proving that dim(Mf t ) > 
\T\D— 1. Since M/ t is a proper face of IPmaghp , we know that dim(Mf t ) < dim(I Pmaghp)- 
So, dim(Mft) = \F\D - 1 and thus, M/ t is a facet of IPmaghpq ■ 

Finally, we consider the set 


N f .ft = {{y,z) G IPmaghp ■ Vjt ~ z f > t = 0}, for some (/',/) G C,t G T. 
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For all / 6 {1, .... |^|} there are t -T f + 1 distinct solutions from the matrices of Figures 
C-l and C-2 which do not belong to N rft . For each of these rows replace the 0’s in the 

columns corresponding to yfr,t < t' < T ft with l’s. 

The remaining matrix will have \T\D affinely independent rows, proving that dim(N rft ) > 
\P\D — 1 . Since Nff t is a proper face of IPmaghp > we know that dim(Nff t ) < 
dim(I P\j agh p) ■ So, dim(N rft ) = \7\D - 1 and thus, N rfl is a facet of IPmaghp- 
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Figure C-l: First Half of Matrix of Solutions to IPmaghp- 
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